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Abstract 

This  thesis  represents  a  methodology  for  the  identi¬ 
fication  of  cost  effective  strategic  force  mixes.  The  meth¬ 
odology  makes  use  of  Response  Surface  Methodology  (RSM) , 
economic  production  functions,  economic  theory,  determin¬ 
istic  models,  and  Lagrangian  techniques  to  identify  cost 
effective  choices.  The  methodology  fits  economic  production 
functions  to  the  response  surface  of  a  nuclear  exchange 
model  (a  linear  programming  problem)  using  RSM.  It  then 
maximizes  these  economic  production  functions  subject  to  a 
cost  constraint  using  the  Lagrangian  technique.  The  use  of 
economic  production  functions  in  this  manner  gives  economic 
insight  into  the  problem  and  results  in  the  development  of 
some  simple  buy  decision  rules  for  determination  of  cost 
effective  force  mixes.  The  classical  use  of  polynomial 
models  does  not  provide  the  same  degree  of  information  as 
the  economic  production  functions,  and  information  gained 
from  the  polynomial  requires  much  more  work. 


A  METHODOLOGY  FOR  IDENTIFYING  COST 


EFFECTIVE  STRATEGIC  FORCE  NIXES 


I. 


A.  OVERVIEW 

There  has  been  a  great  deal  of  effort  placed  on  quanti¬ 
fying  military  weapon  decisions,  the  intent  being  to  maxi¬ 
mize  some  (Measure  of  merit,  usually  destruction.  Often¬ 
times,  cost  is  only  included  as  an  afterthought  in  these 
analyses,  and  not  as  a  key  ingredient  in  making  proper 
optimal  weapon  buys.  In  this  thesis,  a  nuclear  exchange 
model  which  maximizes  the  damage  expectancy  of  a  number  of 
different  types  of  weapons  will  be  analyzed.  The  analysis 
will  make  use  of  economics.  Response  Surface  Methodology 

(RSM),  and  the  properties  of  deterministic  models  to  iden¬ 
tify  cost  effective  strategic  force  mixes.  In  essence,  this 
combination  of  economics,  RSM,  and  deterministic  modelling 
represents  a  methodology  for  identifying  cost  effective 
choices  with  application  to  a  broad  range  of  problems. 

The  thesis  explores  the  problem  of  identifying  cost 
effective  choices  in  a  very  general  manner.  The  intent  is 
that  the  results  can  be  easily  applied  to  other  problems. 
The  specific  application  of  the  methodology  in  this  thesis 
is  to  a  linear  programming  problem  called  a  nuclear  exchange 


Model,  but  any  deterministic  model  could  use  RSM  and 
economics  to  perform  an  analysis  of  cost-effective  choices. 


B.  BACKGROUND 

Response  Surface  Methodology  (RSM)  came  into  being  as  a 
theory  in  the  early  1950’ s.  Through  the  50* s  and  early  60* s 
much  work  was  done  in  developing  RSM  into  a  tool  for  study¬ 
ing  the  results  of  an  experiment.  Through  RSM  an  experimen¬ 
ter  could  fit  a  postulated  surface  to  his  experimental 
results  and  thereby  gain  an  understanding  into  possible 
mathematical  relationships  between  the  factors  in  his  exper — 
iment.  Additionally,  optimal  values  of  the  surface  could  be 
located  using  various  search  techniques. 

While  RSM  has  been  used  extensively  in  experimentation, 
there  has  been  little  application  to  the  results  of  detei — 
ministic  models.  In  this  context,  a  deterministic  model  will 
be  any  MKidel  that  requires  more  than  one  functional  rela¬ 
tionship  to  determine  the  result.  An  example  is  a  Linear 
PrograMuning  problem  that  has  an  objective  function  and  one 
or  Miare  linear  constraints.  Deterministic  models  yield 
results  quite  different  than  experimentation.  For  example, 
there  is  no  sampling  error  in  the  results  of  the  determinis¬ 
tic  model.  That  is,  an  experiment  will  yield  slightly 
different  results  each  time  it  is  repeated,  this  is  called 
saMipling  error,  but  a  deterministic  model  will  yield  the 


same  result  every  tiNte 


The  output  of  a 


it  is  applied 


deterministic  model  is  the  true  surface  in  as  much  as  the 
model  conforms  to  the  reality  it  is  intended  to  simulate. 
Unfortunately,  many  deterministic  models  take  an  incredibly 
long  time  to  calculate,  even  for  a  computer.  The  model  may 
also  represent  so  many  combinations  of  factors,  that  it  is 
not  worthMhile  to  calculate  the  true  value  for  each.  RSM 
represents  a  method  to  "see"  the  resulting  surface  with  a 
minimal  of  model  runs. 

Since  the  deterministic  model  is  itself  the  true  model 
of  the  phenomena,  any  postulated  mathematical  equation  (also 
called  a  model)  will  not  fit  the  surface  perfectly  and  in 
effect  will  be  biased.  Bias  represents  systematic  error  and 
is  error  caused  by  an  incorrect  model  being  fitted  to  the 
data.  A  major  defficiency  currently  in  RSM  is  the  lack  of 
experimental  designs  that  minimize  this  bias.  An  experi¬ 
mental  design  is  an  "optimal"  set  of  factor  values  which 
minimize  the  number  of  observations  required  to  minimize  the 
variance  error  in  fitting  a  functional  form  to  the  results. 
For  instance,  if  one  selects  (postulates)  a  straight  line  as 
his  model  that  would  best  fit  his  results  (response),  only 
two  observations  are  needed  to  fit  a  line  to  the  response. 
This  represents  the  smallest  design  needed  to  fit  the  model. 
The  theory  of  experimental  design  allows  one  to  select  the 
least  number  of  points  needed  to  accurately  fit  the  model. 
Most  of  the  designs  used  in  RSM  seek  to  minimize  variance 
error  and  do  not  specifically  take  bias  into  account. 
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1.3 


In  fitting  a  model  to  a  surface.  Ordinary  Least  Squares 
(OLS)  is  the  most  common  approach.  The  Sue  of  Squares  of 
Error  (SSE)  is  the  total  squared  error  between  the  postu¬ 
lated  model  and  the  actual  surface.  SSE  is  a  combination 
of  variance  error  and  bias  error.  It  is  possible  through 
the  use  of  replications  to  find  out  what  portion  of  the 
error  is  due  to  lack  of  fit  (bias)  or  sampling  (variance). 
In  the  application  of  models  to  the  surface  of  a  determinis¬ 
tic  model,  the  SSE  represents  all  bias  error.  Therefore, 
the  model  that  most  closely  represents  the  true  surface  will 
be  that  which  minimizes  SSE. 

Most  mathematical  models  postulated  to  a  reponse  sur — 
face  are  first  and  second  order  models.  While  models  of 
this  form  often  give  very  good  “fit'*  to  the  responses,  the 
relationship  among  the  coefficients  of  the  various  factors 
can  be  obscure  due  to  magnitude  differences  in  the  inputs  or 
difficulty  in  interpreting  the  meaning  between  linear  coef¬ 
ficients,  quadratic  coefficients,  and  interaction  coef¬ 
ficients  of  the  same  factor.  The  reason  a  second  order 
polynomial  gives  such  a  good  fit  is  that  it  represents  "an 
approximating  second  order  Taylor  expansion  of  the  true  but 
unknown,  response  surface”  (Box  and  Draper,  1971:734).  An 
example  of  a  second  order  polynomial  is 

Y  -  2  «•  6X  +  3Y  -  2XY  +  3X2  +  Y2. 

There  are  other  models  which  may  lend  themselves  better  for 
interpretation  of  the  coefficients,  but  how  good  would  they 


graphically  showing  the  effect  of  interactions  on  the  re¬ 
sponse.  Otto  and  Werner,  1983,  used  three-dimensional  plots 
to  show  two-way  interactions  effect  on  certain  chemical 
reactions  (Otto  and  Werner,  1983:246,  263).  Frazer  and 
Brand,  1980,  demonstrated  the  use  of  graphical  techniques 
for  data  analysis  using  RSM. 

Frazer  and  Brands’  paper  is  an  important  divergence 
from  the  classic  chemical  analysis  using  RSM.  Frazer 
attempts  to  fit  the  the  response  surface  to  a  computer  model 
of  the  surface.  The  abstract  explains  the  intent  of  his 
paper.  Frazer  uses 

...  a  mathematical  nonlinear  least-squares  tech¬ 
nique  that  fits  the  entire  experimental  data  set 
to  the  different  kinetic  models.  We  also  show 
that  when  these  NLSS  mathematical  techniques  are 
used  directly  to  test  kinetic  models,  estimates  of 
physical  constants  are  improved,  fewer  graphical 
representations  are  required,  and  a  most  likely 
mechanism  is  more  easily  and  dependably  identi¬ 
fied"  (Frazer  and  Brandt,  1983:  1730). 

Frazer  is  comparing  the  response  surface  of  actual  experi¬ 
mental  data  with  the  response  surface  generated  by  nonlinear 
models.  Using  these  nonlinear  deterministic  models,  he  is 
able  to  more  closely  approximate  the  real  underlying  rela¬ 
tionship.  While  there  are  some  similarities  to  the  thesis 
of  this  paper,  Frazer’s  goal  is  quite  different;  he  is 
trying  to  estimate  an  equation  to  a  set  of  experimental  data 
that  contains  measurement  error. 

In  1964,  Wu  authored  a  series  of  articles  using  RSM  to 
predict  tool  wear.  Whether  it  was  a  result  of  these  arti- 


analysis  o-f  variance  test  (ANOVA)  when  the  model  is  simple, 
and  the  modified  principal  component  analysis  otherwise 
(Snee,  1972:61).  The  only  article  regarding  analysis  of 
response  curves  for  deterministic  models  is  Smith,  1979. 
Smith  used  RSM  to  perform  a  multi -dimensional  parametric 
analysis  of  a  model  rather  than  the  classical  one  dimen¬ 
sional  sensitivity  analysis. 

C.  APPLICATION  OF  RESPONSE  SURFACE  METHODOLOGY 

RSM  has  been  applied  to  numerous  fields  of  science. 
Almost  any  three  months  of  publication  would  yield  a  dozen 
applications.  In  this  section  three  main  applications  will 
be  discussed,  the  first  as  an  example  of  the  historical 
development  of  RSM,  and  the  other  two  as  they  illustrate  the 
broad  range  of  RSM  and  potentially  apply  to  this  thesis  sub¬ 
ject. 

RSM  has  found  a  great  deal  of  acceptance  in  the  chemi¬ 
cal  industry.  Historically,  this  is  where  the  first  appli¬ 
cations  of  RSM  were  employed.  Today,  there  has  been  no 
abatement  in  the  use  of  RSM  in  analyzing  chemical  experi¬ 
ments.  A  recent  application  is  Otto,  Rentsch,  and  Werner, 
1983  M.  .  .  optimal  conditions  for  determination  of  CO<II) 
are  deduced  from  response  surface  studies,  considering  the 
sensitivity  and  the  blank  absorbance  as  responses"  (Otto  and 
others,  1983:267).  The  chemical  industry  has  used  RSM  for 
more  than  predicting  responses,  but  also  as  a  means  of 


be  used  to  simulate  a  reel  life  process,  end  provided  the 
computer  eodel  does  not  employ  any  monte-carlo  techniques, 
the  result  Mill  always  be  the  same.  In  other  mords,  the 
computer  has  no  measurement  errors  associated  with  its  ex¬ 
periment,  which  means  no  variance.  This  requires  the  use  of 
minimum  bias  designs.  Very  little  work  has  been  done  in 
this  area  of  applying  RSM  to  deterministic  models.  Smith, 
1979,  is  one  exception.  The  payoff  of  applying  RSM  to 
deterministic  models  is  great  considering  the  tremendous  run 
time  of  some  computer  models.  To  be  able  to  predict  the 
computer  result  through  the  use  of  some  mathematical  equa¬ 
tion  also  provides  insight  into  the  problem. 

While  RSM  is  an  important  tool  in  the  estimating  of  re¬ 
sponses,  to  stop  here  would  show  a  lack  of  appreciation  for 
the  full  theory. 

It  is  of  utmost  importance  that  the  experimenter 
be  able  to  appreciate  the  results  of  an  RSM  study. 

This  point  cannot  be  overemphasized.  There  have 
been  cases  where  rather  elaborate  mathematical 
models  have  been  proposed,  the  appropriate  data 
collected,  the  models  fitted,  and  the  results 
presented  in  the  form  of  mathematical  equations, 
only  to  have  the  practical  implications  of  these 
results  remain  unappreciated  and  therefore  unex¬ 
ploited  (Hill  and  Hunter,  1966(573). 

There  are  a  number  of  good  articles  on  the  analysis  of 
response  surfaces.  The  three  most  referenced  articles  are 
Box,  1954;  Box  and  Youle,  1955;  and  Snee,  1972.  Snee  re¬ 
views  a  couple  of  practical  applications  using  principal 
component  analysis,  a  modified  principal  component  analysis, 
and  a  shape  curve  analysis.  Snee  recommends  a  general 


2. 


as  in  controlling  the  environment  of  the  experiment,  instru¬ 
ment  errors,  etc.),  while  the  bias  error  measures  error 
caused  by  specifying  an  incorrect  model  to  tha  surface.  In 
most  experimental  design  work,  the  focus  has  been  on  mini¬ 
mizing  the  variance  error  with  little  done  to  minimize  the 
bias  error.  "Initially,  criteria  for  judging  goodness  of 
designs  were  largely  concerned  with  variance,  ....  The 
question  of  bias  due  to  inadequacy  of  the  ...  [model  3  was 
given  somewhat  secondary  consideration"  (Karson  and  others, 
1969:461,  462).  Very  little  work  has  been  done  in  this  area 
of  developing  designs  to  minimize  bias  error.  Karson, 
Manson,  and  Haders’  work  in  1969,  Karson,  1970,  and  Box  and 
Draper,  1963,  are  the  only  articles  located  on  this  subject. 
Myers  summarizes  Box  and  Draper’s  work  in  one  chapter  in  his 
book  (Myers,  1976:196-218). 

The  reason  that  emphasis  has  been  placed  mainly  on 
minimizing  error  due  to  variance  is  that  most  of  the  appli¬ 
cations  of  RSM  have  been  experiments.  The  most  common  early 
experiments  were  performed  using  different  chemicals  and 
environments  to  test  reaction.  In  the  development  of  RSM, 
it  was  very  important  to  minimize  the  error  due  to  measure¬ 
ment  of  the  response  in  the  laboratory. 

Being  able  to  minimize  bias  becomes  a  crucial  issue 
when  one  moves  from  the  application  of  RSM  to  experiments 
using  imperfect  measuring  instruments  to  the  application  of 
computer  generated  responses  to  problems.  The  computer  can 
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FIGURE  2.3 


In  analyzing  the  coef -ficients  of  a  hypothesized  model, 
it  is  very  useful  to  have  coefficients  which  are  uncorre¬ 
lated;  in  that  way  the  error  is  not  confounded.  Confounded 
error  eeans  that  it  can  not  be  distinguished  which  factors 
were  the  sain  contributors  of  the  error,  such  as  an  interac¬ 
tion  or  quadratic  term.  In  I960  Box  and  Behnken  suggested 
som  possible  new  designs,  "All  the  designs  we  discuss 
possess  a  high  degree  of  othogonality;  in  fact,  only  the 
constant  term  .  .  .  and  the  quadratic  estimates  ...  are 
correlated  one  with  another"  (Box  and  Bhenken,  1960*457). 
These  designs,  developed  by  Box  and  Bhenken,  are  employed  in 
this  study  and  are  referred  to  as  variance  designs;  that  is, 
designs  which  seek  to  minimize  variance. 

The  total  error  in  a  model  can  be  divided  into  error 
due  to  variance  error  and  error  due  to  bias  error.  The 
variance  error  reflects  measurement  error  of  some  kind  (such 


bined,  enable  the  ixpcriiMntar  to  Mki  an 
efficient  empirical  exploration  of  the  system  in 
which  ho  is  intsrsstsd  (Myers,  1976*1). 

Box  and  Wilson  first  dsvslopsd  ths  thsory  of  rssponss 
surfacs  aethodolgy  in  19S1.  Ths  thsory  continusd  to  bs 
sxpandsd  through  ths  50’ s  and  oarly  60’ s  Mith  ths  devel- 
opssnt  of  new  experimental  dssigns  (Box  and  Hunter,  19S7| 
Box  and  Draper,  1963).  There  are  a  nusbsr  of  different 
types  of  dssigns  with  varying  degress  of  usefulness,  depend¬ 
ing  on  the  type  eodel  being  fitted.  The  eost  common  eodels 
are  first-order  and  second-order .  A  first-order  eodel  is  a 
eodel  of  linear  fore)  a  second-order  eodel  adds  quadratic 
variables  as  well  as  two-way  interactions.  A  first  order 
eodel  can  also  have  interaction  teres.  The  following  equa¬ 
tions  illustrate  these  models. 

Y  ■  pQ  +  PjXj  +  P2X2  (First-order  eodel) 

Y  -  “o  ♦  '1X1  *  *2X2  +  '.2V2  *  »nx?  *  »22XI 

(Second-order  eodel) 

A  design  for  a  first-order  eodel  is  called  a  first-order 
design  and  likewise,  a  second-order  design  for  a  second- 
order  eodel.  Myers,  1976,  gives  a  good  overview  of  the 
different  type  of  designs. 

An  exaeple  of  an  experieental  design  is  illustrated  in 
figure  2.3.  The  einieum  and  maximum  levels  to  be  tested  of 
the  cheeical  are  represented  by  MIN  and  MAX.  The  design  is 
a  simple  uniform  design  measuring  three  different  levels  of 
each  chemical  for  a  total  of  nine  combinations. 


fans  that  ths  Modal  is  biassd 


It  Mill  not  -fit  ths  data 


perfectly  due  to  this  bias,  the  inability  to  completely 
specify  the  factors  involved  in  the  eodel.  Variance  error, 
the  other  type  of  error  in  an  experiment,  is  caused  by  the 
inability  to  hold  certain  factors  constant  in  the  experiment 
and  to  perfectly  measure  the  response.  For  both  of  these 
reasons,  the  model  Mill  usually  not  fit  the  responses  as 
sell  as  in  figure  2.2.  Using  a  technique  such  as  least 
squares  alloMs  the  determination  of  the  error  that  is  due  to 
poor  measurement  of  the  responses  or  error  due  to  an  in¬ 
correct  model  postulated  by  the  experimenter. 

RSM  helps  in  the  design  of  the  experiment  in  order  to 
determine  the  least  number  of  measurements  needed  to  fit  the 
model  to  the  data.  RSM  defines  the  statistics  used  in 
calculating  horn  accurate  the  postulated  model  is  to  the  real 
data.  If  a  mathematical  model  can  be  estimated  using  RSM  to 
explain  the  process  Mith  a  high  degree  of  accuracy,  then  one 
no  longer  needs  to  perform  the  experiment  to  find  out  the 
response;  responses  are  calculated  using  the  mathematical 
model  (for  the  combinations  of  factors  Mithin  the  limits  of 
the  experiment).  Lastly,  RSM  uses  techniques  such  as  steep¬ 
est  ascent  to  find  the  combinations  for  optimal  output. 

Response  8urface  Methodology  .  .  .  is  essentially 
a  particular  set  of  mathematical  and  statistical 
methods  used  by  researchers  to  aid  in  the  solution 
of  certain  types  of  problems  Mhich  are  pertinent 
to  scientific  or  engineering  processes.  ...  The 
response  surface  procedures  are  a  collection  in¬ 
volving  experimental  strategy,  mathematical  meth¬ 
ods,  and  statistical  inference  Mhich,  when  com- 


FIGURE  2.1 


FIGURE  2.2 


tried,  and  the  response  (number  of  bushels)  is  Measured  for 
each  combination.  Next,  the  scientist  Mill  attenpt  to  fit  a 
postulated  mathematical  model  to  the  responses. 

Figure  2.1  is  a  graphic  picture  of  the  example.  Four 
different  combinations  of  fertilizer  amount  and  depth  of 
planting  have  been  employed.  The  response  for  each  has  been 
measured  and  plotted  on  the  graph  (the  z  axis  represents 
yield,  the  x  and  y  axis,  fertilizer  amount  and  depth).  The 
responses,  marked  Mith  an  asterisk,  represent  the  response 
surface.  Using  the  method  of  least  squares,  a  math  model 
can  be  fitted  to  the  response  surface  as  in  figure  2.2. 

Obviously,  there  are  many  more  factors  involved  in  the 
yield  of  a  crop  then  just  fertilizer  and  depth  of  the  cut. 
Other  factors  mould  be  the  amount  of  rain,  amount  of  sun, 
type  of  soil,  plant  proximity  to  each  other,  weed  influence, 
etc.  The  exclusion  of  these  other  factors  in  the  model 


A.  INTRODUCTION 


There  is  sxtsnsivs  literature  on  Response  Surface 
Methodology  (RSM) ,  dstarministic  nodal s,  and  production 
functions  in  economics,  but  littls  in  applying  all  thrss 
together.  Ths  topic  of  this  study,  a  methodology  for  iden¬ 
tifying  cost  effective  strategic  force  sixes,  uses  a  combi¬ 
nation  of  disciplines  apparently  unique  to  the  scientific 
literature.  Still,  this  combination  of  RSM,  economics,  and 
deterministic  models  suggests  four  potential  areas  in  Mhich 
to  review  the  literatures  Response  Surface  Methodology, 
application  of  RSM  to  deterministic  models,  economic  appli¬ 
cations  using  RSM,  and  economic  production  models. 

B.  THEORY  OF  RESPONSE  SURFACE  METHODOLGV 

Response  Surface  Methodology  is  a  method  for  fitting 
mathematical  models  to  surfaces  generated  by  experiment. 
These  mathematical  models  are  equations  that  the  scientist 
has  postulated  will  describe  the  response  produced  by  his 
experiment.  For  example,  suppose  one  wishes  to  know  the 
yield  of  a  particular  crop  depending  on  the  amount  of  ferti¬ 
lizer  and  the  depth  of  the  seed  planted.  The  response  of 
the  experiment,  the  yield  of  the  crop,  can  be  measured  in 
number  of  bushels.  Than,  a  number  of  different  combina¬ 
tions  of  amount  of  ferilizer  and  depths  of  seed  planted  are 


lowd  by  a  thrac  factor  aodal.  Tha  fol lotting  functions  ars 
aatiiutad  for  ths  surfacst  1)  Polynomial,  2)  Cobb-Dougl as, 
3)  Generalized  CES,  4)  Multiplicative  VES,  5)  Additive  VES. 
Ths  rssults  of  sach  aodal  ars  coaparsd  and  ths  fit  of  ths 
aodel  tsstad.  Also,  rslationships  in  ths  coafficiants  ars 
SMaainadi  and  ths  aiodsls  svaulatsd  to  sss  what  kind  of 
intsrpratation  can  bo  givsn  to  sach  aodal. 

Finally,  a  cost  constraint  is  addsd  to  ths  deter- 
ainistic  aodel.  Again,  thsrs  ars  two  approachas.  Ons,  ths 
cost  constraint  can  ba  addsd  within  ths  datarainistic  aodal. 
Two,  it  can  bo  ussd  out si da  of  ths  aodal.  Ths  sscond 
approach  would  aaxiaizs  ths  ostiaatad  function  subject  to  a 
cost  constraint  using  ths  Lagrangian  asthod.  Ths  goal  is  to 
idsntify  cost  offsctivo  choicss. 

In  suaaary,  two  aajor  goals  ars  to  analyzs  ths  rela- 
tionships  aaong  ths  diffsrant  functions  and  to  sse  how 
ussful  ths  additional  sconoaic  intsrpratation  is  in  giving 
the  decision  aaker  mors  information,  specifically,  informa¬ 
tion  that  helps  in  the  deteraination  of  cost  affective 
aixes. 
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tween  the  two  typos  of  oodals  and  aconooic  intorprotation  of 
tha  results.  Finally,  add  cost  to  those  functions  to  sea 
what  typo  of  cost  effective  results  are  given. 

F.  APPROACH 

As  a  scans  of  answering  the  questions  that  have  been 
raised  in  this  proposal,  a  deterministic  model,  with  mili¬ 
tary  aplications,  will  be  used  to  test  these  ideas.  The 
deterministic  model  is  a  nuclear  exchange  model  that  com¬ 
putes  the  optimal  damage  for  a  given  number  of  weapons  and 
weapon  types.  This  is  a  common  production  problem  with 
output  measured  in  damage  expectancy  and  the  input  variables 
are  a  given  number  of  weapons  per  weapon  type. 

The  first  step  is  to  research  the  nuclear  exchange 
model  to  fully  understand  what  exactly  it  is  doing.  What 
type  of  model  is  it?  Does  it  contain  any  critical  points? 

Next,  research  is  accomplished  to  find  appropriate 
experimental  designs  with  an  emphasis  on  those  which  mini¬ 
mize  bias.  There  is  very  little  work  done  in  this  area,  but 
it  may  be  possible  to  use  some  present  design  which  may  have 
bias  minimizing  potential.  Box  and  Draper,  1959,  is  a 
potential  resource  in  this  area.  Additionally,  a  nonlinear 
estimation  technique  is  required  to  fit  the  CES  production 
function  to  the  data.  All  other  models  are  estimated  using 
ordinary  least  squares. 

After  deciding  on  a  design  for  each  function,  the 
computer  runs  are  made.  A  two  factor  model  is  used,  fol- 
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D.  RESEARCH  QUESTION 


Can  econoaic  interpretations  ba  obtained  froa  the  typi¬ 
cal  linear  deterainistic  aodel  and  its  response  surface 
aodel?  How  does  this  differ  froa  coaaon  experiaental  pro¬ 
duction  functions?  Can  RSM  and  experiaental  design  be  used 
to  fit  production  functions,  such  as  the  Cobb— Douglas,  a 
generalized  CES,  a  Multiplicative  VES,  and  an  additive  VES, 
to  a  response  surface?  How  are  these  aodels  interpreted 
and  what  relationships  exist  between  the  various  aodels, 
including  the  second  order  polynoaials?  Which  aodels  fit 
the  response  surface  best,  and  why  is  this  so?  What  is  the 
tradeoff  between  a  functional  fora  that  predicts  well  and 
one  that  is  easy  to  interpret  and  analyze?  What  happens 
when  cost  is  added  to  the  aodel,  what  type  of  answers  will 
the  response  surface  give?  Does  the  use  of  econoaic 
production  functions  provide  greater  insight  into 
deteraining  cost  effective  aixes?  How  can  others  use  this 
application  in  their  analyses? 

E.  OBJECTIVES  QF  IHg  RESEARCH 

The  objectives  of  the  research  are  to  use  RSM  and 
appropriate  experiaental  designs  to  apply  coaaon  econoaic 
production  functions  to  a  response  surface.  Then  evaluate 
these  production  functions  with  second  order  equations  over 
a  range  of  factor  levels  (nuaber  of  variables  in  the  equa¬ 
tion).  Coapare  the  results  and  look  for  relationships  be- 


being  introduced  into  the  problem.  If  the  SSE,  which  is  all 
bias  error,  can  be  reduced  to  a  sufficiently  seall  per¬ 
centage,  then  insight  into  the  problen  can  be  gained  by 
studying  the  functional  fore,  at  to  mention  the  savings 
from  being  able  to  predict  results  without  having  to  rerun 
the  deterministic  model. 

Finally,  the  most  important  interpretation  of  functions 
using  RSil  will  be  those  involving  cost.  How  are  the  results 
of  the  model  modified  by  the  addition  of  a  cost  restraint? 
There  are  two  approaches  in  adding  cost  to  the  analysis.  A 
cost  restraint  could  be  included  in  the  deterministic  model, 
or  the  estimated  production  functions  could  be  maximized 
subject  to  a  cost  restraint.  The  question  is  then,  where 
does  one  include  cost  in  the  analysis,  i.e. ,  as  a  restraint 
or  as  another  function  to  sub ject  the  model  to? 

C.  PROBLEM  STATEMENT 

It  appears  that  little  has  been  done  to  evaluate  in 
economic  (or  other  nonstatistical )  terms  the  models  used  to 
capture  the  response  surface  of  a  deterministic  model.  What 
economic  insights  are  revealed  by  response  surface  modelling 
and  how  do  these  compare  with  results  of  experimental  models 
of  similar  economic  phenomena?  What  information  can  be 
gained  for  identifying  cost  effective  choices?  A  methodol¬ 
ogy  is  needed  that  evaluates  and  integrates  RSM,  determin¬ 
istic  modeling,  and  economics  and  into  a  form  that  will  be 
useful  to  others  using  RSM  in  analysis. 
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fit  the  rnponM?  Also  which  is  sots  important  in  s  re¬ 
sponse  surface  equation,  predicting  ability  or  interpreting 
ability?  Obviously  these  answers  would  depend  on  the  appli¬ 
cation,  but  it  appears  to  be  an  area  that  has  had  little 
research. 

There  are  other  functions,  such  as  those  used  in  eco- 
nonics,  which  could  also  be  postulated  as  eodels.  They  are 
the  Cobb-Oouglas,  the  generalized  Constant  Elasticity  of 
Substitution  (CES) ,  a  Variable  Elasticity  of  Substitution 
(VES)  Multiplicative  model,  and  a  VES  additive  model. 

This  last  set  of  models,  which  have  been  used  in  eco¬ 
nomics,  may  have  some  advantage  over  a  second  order  poly¬ 
nomial  in  that  the  coefficients  have  easily  computed  eco¬ 
nomic  interpretations,  such  as  marginal  products,  elastici¬ 
ties,  etc.  Other  advantages,  at  least  in  the  case  of  the 
Cobb-Douglas  and  CES  production  function  are  that  they  are 
not  affected  by  magnitude  differences  among  the  input  vari¬ 
ables.  The  VES  models  contain  interaction  terms  like  a 
second  order  polynomial,  and  the  functional  expressions 
needed  to  calculate  basic  economic  values  are  much  simpler 
than  the  polynomial.  Whether  these  models  can  predict  as 
well  as  a  second  order  polynomial  is  an  interesting  ques¬ 
tion. 

No  matter  what  the  deterministic  model  is,  the  true 
surface  very  rarely  will  ever  be  captured  perfectly  by  a 
mathematical  equation.  From  the  very  start,  bias  error  is 


cles  or  not,  the  manufacturing  industry  has  made  Hide  appli¬ 
cation  to  various  tool  near  problems.  Wu  applies  a  couple 
of  different  models  to  approximate  tool  near  based  on  five 
input  factors:  speed,  feed,  depth  of  cut,  angle,  and  nose 
radius.  His  articles  are  of  special  interest  because  he 
uses  a  multiplicative  model  to  estimate  tool  life.  This 
model  is  very  similar  in  form  to  a  Cobb-Douglas  (C-D)  pro¬ 
duction  function,  which  is  used  in  economics  extensively  and 
was  used  in  this  thesis.  The  Cobb-Douglas  production  func¬ 
tion  provides  economic  insight  and  interpretation  to  a  proc¬ 
ess  which  would  be  much  more  difficult  using  an  ordinary 
second  order  model . 

Wu’s  use  of  a  multiplicative  function  like  the  C-D  was 
not  to  gain  any  economic  insight,  but  rather  because  pre¬ 
vious  investigators  had  proposed  a  multiplicative  functional 
relationship  (Wu,  1964:106).  Wu  used  RSI1  to  estimate  a 
multiplicative  function  for  the  purpose  of  saving  time  and 
money  being  spent  on  tool  life  measurement. 

A  third  application  of  RSM  is  in  nutrition.  Two  arti¬ 
cles  which  discuss  the  use  of  RSM  in  feeding  chicks  are 
Toyomizu,  Akiba,  Horiguchi,  and  Matsumoto,  1982,  and  Roush, 
Petersen,  and  Arscott,  1979.  These  serve  as  excellent  ex¬ 
amples  of  the  varied  use  of  RSM,  Toyomizu,  et.  al . ,  used  a 
creative  triangular  design  to  estimate  a  second  order  equa¬ 
tion  with  three  variables  (Toyomizu  and  others,  1982:887). 
Both  of  these  articles  applied  RSH  to  determining  the 


opt i Ml  combination  of  inputs  to  maximize  the  output. 

There  have  bem  many  other  applications  of  RSM.  These 
three  were  mentioned  because  they  give  an  idea  of  the  broad 
range  of  application  of  RSM. 

D.  ECONOMIC  ANALYSIS  USING  RESPONSE  SURFACE  METHODOLOGY 

The  specific  use  of  RSM  in  economics  is  scarce. 
Because  of  the  nature  of  economics,  specific  applications 
would  have  to  focus  on  developing  large  scale  economic 
models,  either  deterministic  or  simulation,  and  then  fit 
economic  functions  to  these  responses.  Only  in  these  appli¬ 
cations  would  economics  be  fully  utilizing  RSM. 

Burdick  and  Naylor  discuss  a  potential  application 
using  Klein’s  six-equation  econometric  model  of  the  United 
States  and  a  Cobb -Doug las  utility  function  (Burdick  and 
Naylor,  1969:29-31).  Unfortunately,  they  do  not  actually 
apply  RSM{  they  only  discuss  how  it  could  be  used.  This  was 
the  only  article  found  on  a  specific  deterministic 
application  of  RSM  in  economics. 

Karlinger  and  Attanasi ,  1980,  use  RSM  to  evaluate  the 
effect  of  the  risk  of  a  flood  on  a  person’s  willingness  to 
buy  flood  insurance.  They  used  a  simulation  program  to 
model  their  problem  and  then,  fit  a  utility  function  to  the 
surface.  This  was  the  only  other  specific  application  of 
RSM  in  economics  found. 

The  most  likely  reason  for  the  lack  of  attention  in 
economics  to  this  area  is  that  a  model  (deterministic  or 
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simulation),  is  only  an  attempt  to  capture  the  real  world, 

and  the  -function  -fitted  to  the  response  surface  of  the  model 

is  an  attempt  to  capture  this  surface.  The  potential  for 

error  is  doubled  with  the  use  of  two  models. 

Now  there  are  two  response  surfaces  we  can  con¬ 
sider:  one  in  the  real  world  relating  the  real 

world  response  to  real  world  factors  and  the  other 
is  the  corresponding  response  surface  in  the 
model.  It  is  our  hope  that  the  two  will  be  very 
close,  but  this  is  much  more  likely  to  be  true 
near  the  real  world  data  points  than  far  away  from 
them.  If  our  search  leads  us  to  an  optimum  point 
far  from  the  real  world  data  points,  any  resulting 
conclusions  about  the  location  of  the  optimum 
point  in  the  real  world  should  be  made  with  great 
caution  (Burdick  and  Naylor,  1969:31). 

While  economics  has  not  made  many  specific  applications 
of  RSN,  it  has  sought  to  capture  data  (real  world  data)  in 
functional  forms.  Since  the  data  are  real  world,  it  does 
not  come  in  the  experimental  design  form  used  in  RSM.  It 
is  not  experimental,  but  rather  observation  related. 
Similar  to  RSN,  the  techniques  of  least  squares  is  widely 
used  to  fit  models  to  real  world  data.  These  real  world 
data  are  in  a  sense  a  response  surface.  The  only  difference 
is  that  it  is  random  and  not  generated  from  an  experimental 
design. 

Therefore,  a  second  approach  in  performing  a  literature 
review  was  to  locate  articles  about  well  known  production 
functions  used  in  economics.  This  will  form  a  nucleus  of 
information  that  can  be  referenced  later  after  attempting  to 
fit  these  models  to  deterministic  data. 


Thera  is  a  lot  o-f  literature  concerning  production 
functions  in  economics.  Those  with  a  good  overall  descrip¬ 
tion  include  Ferguson,  1969;  Henderson  and  Quandt,  1980; 
Nicholson,  1978;  and  Silberberg,  1978. 

Zellner,  1966,  deals  specifically  with  the  Cobb-Douglas 
(C-D)  production  function.  The  general  fore  of  the  C-D  for 
two  variables  is  given  in  equation  2.1. 

Y  ■  pK*L*—1  <2. 1) 
For  the  purpose  of  explanation,  consider  the  variables  K  and 
L  to  represent  capital  and  labor.  The  powers  of  the  var¬ 
iables  K  and  L  are  the  output  elasticities  of  capital  and 
labor.  These  measure  the  percent  change  in  output  (produc¬ 
tion)  for  a  one  percent  change  in  an  input  (either  capital 
or  labor).  As  an  example,  if  a  equaled  0.5,  then  a  one 
percent  increase  in  capital  would  cause  a  0.5  percent 
increase  in  production.  The  a  and  l-«  terms  are  also  a 
measure  of  the  input  intensities.  The  greater  the  a  value, 
the  more  capital  intensive  the  production  process  is  (Brown, 
1968x48). 

The  "p"  value  in  equation  2.1  measures  technical  prog¬ 
ress.  Technical  progress  is  an  abstract  concept  capturing 
the  importance  of  technology  in  the  production  process. 

Arrow,  Chenery,  Minhas,  and  Solow  first  introduced  a 
more  general  form  of  the  C-D,  a  constant  elasticity  of 
substitution  (CES)  function,  in  1961.  Uzawa,  1962,  gener¬ 
alized  it  for  more  than  two  parameters.  Constant  elasticity 


of  substitution  means  that  the  ability  to  substitute  among 
the  inputs  and  maintain  constant  output  does  not  change  -for 
varying  levels  of  the  inputs.  Equation  2.2  is  an  example  of 
the  CES  for  two  inputs. 

— P  —p  —  V/P 

Y  -  pC«K  +  (1-«)L  3  (2.2) 

The  variable  p  is  an  efficiency  (or  technology)  parame¬ 
ter}  «  is  a  distribution  parameter}  P  is  a  substitution 
parameter}  and  V  measures  returns  to  scale. 

The  technology  parameter  has  the  same  abstract  meaning 
as  in  the  Cobb-Dougals  above)  it  is  a  measure  of  the  imp or — 
tance  of  technology  in  the  production  process. 

The  distribution  parameter  measures  the  intensity  of 
each  input  in  the  technology.  If  «  is  large  the  technology 
is  K  intensive.  In  other  words,  the  technology  relies  more 
on  the  input  K  than  L  to  achieve  its  efficiency. 

P  measures  the  ease  of  substitution  among  the  input 
variables  (K  and  L) .  The  ease  of  substitution  is  a  diffi¬ 
cult  parameter  to  conceptualize.  The  higher  the  ease  of 

substitution,  the  easier  it  is  to  substitute  inputs  and 
still  maintain  production,  the  lower  the  ease  of  substitu¬ 
tion,  the  harder  it  is  to  substitute  inputs  for  each  other 
and  maintain  output.  For  high  P  values,  which  implies  that 
the  ease  of  substitution  is  difficult,  the  production  pro¬ 
cess  is  referred  to  as  proportional ,  that  is,  proportionate 
levels  of  each  input  are  needed  for  production.  For 

example  one  cannot  trade  labor  for  capital  very  easily  in 


this  process.  As  a  reference,  the  C-D’s  elasticity  of 
substition  is  one,  which  is  somewhere  between  low  and  high 
ease  of  substitution.  Brown,  1968,  is  an  excellent  source 
for  a  detailed  discussion  of  the  CES  as  well  as  McFadden, 
1963. 

Both  the  C-D  and  CES  are  homogeneous  production  func¬ 
tions.  A  homogeneous  production  function  is  one  where  the 
relationship  given  in  equation  2.3  holds;  such  a  function  is 
homogeneous  of  degree  k. 

f  <tK,  tl_>  *  tkf(K,  L)  (2.3) 

The  most  common  degree  of  homogeneity  is  k  »  1.  In 
less  technical  terms,  a  doubling  of  both  inputs  would  cause 
output  to  double.  If  a  function  is  not  homogeneous,  a 
doubling  of  the  inputs  may  not  even  have  the  same  effect  at 
different  levels  of  the  inputs.  Two  recent  functions  used 
in  economics  which  are  not  homogeneous  are  the  multiplica¬ 
tive  nonhomogeneous  <MNH)  (Vinod,  1972*531-543)  and  the 
additive  nonhomogeneous  (ANH)  (Sudit,  1973*499-514).  Both 
functions  are  an  attempt  to  further  unrestrict  economic 
problems  for  cases  where  variable  elasticities  of  substitu¬ 
tion  (VES)  are  more  appropriate.  Vinod  and  Sudit  applied 
their  production  functions  to  telecommunication  problems 
where  the  ability  to  substitute  among  capital  and  labor 
changes  over  time. 

Equation  2.4  is  the  MNH  function.  The  'A'  is  also  a 
technology  measure  as  in  the  C-D  and  CES.  The  elasticities 


of  output  require  calculation.  The  advantage  of  the  MNH  as 
a  non— hoaoqanaous  production  function  is  that  iBost  of  tha 


aconoaic  values  can  ba  easily  calculated. 

Y  -  AK*+wLN(L>Lp  (2.4) 

Tha  ANH  does  not  have  coefficients  with  siaple 
interpretations,  but  like  the  MNH,  the  output  elasticities, 
the  Marginal  products,  etc.,  can  be  easily  cooputed. 

Later  in  the  thesis,  the  MNH  and  ANH  will  be  explored 
eore  fully  in  a  real  exaeple  of  the  results. 

The  only  article  found  relating  economics  and  response 
surfaces  is  Burdick  and  Naylor,  1969.  Their  article  is  a 
good  article  for  the  economist  attempting  to  understand  the 
theory  of  response  surface  methodology.  It  does  not  give 
any  specific  examples  of  the  application  of  RSM  to 
economics,  although  it  does  suggest  some  applications. 

E.  CONCLUSION 

Response  Surface  Methodology  has  been  around  for  over 
thirty  years  and  has  had  wide  application  in  many  fields. 
Most  applications  have  dealt  with  experiments  having  some 
sampling  error  in  the  data.  Very  little  literature  has  been 
written  about  the  application  of  RSM  to  deterministic 
models.  Detereinistic  models  do  not  contain  error  due  to 
sampling  and  therefore,  require  statistical  designs  that 
minimize  bias  error  as  apposed  to  error  due  to  variance. 
Again,  the  literature  on  minimum  bias  designs  is  slight. 
This  is  to  be  expected  given  the  history  of  application  of 
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RSM  to  industry  experiments.  Additionally,  littls  work  has 
been  dons  in  the  economic  interpretation  of  RSM  models. 

The  economic  interpretation  o-f  RSM  equations  to  deter-* 
mini stic  models  provides  excellent  further  study.  This  work 
would  develop  a  minimum  bias  design  using  the  theory  devel¬ 
oped  by  Karson,  the  fitting  of  common  economic  production 
functions  to  deterministic  data,  and  the  comparison  of  the 
different  models.  In  light  of  this  literature  review,  this 
would  represent  the  tying  together  of  a  number  of  theories 
to  produce  a  useful  document  for  the  study  of  deterministic 
model  responses. 
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A.  INTRODUCTION 


The  methodology  used  to  idanti-fy  cost  eff active  stra¬ 
tegic  -force  sixes  is  composed  of  six  major  steps.  The  steps 
are  as  follows* 

1)  Select  the  deterministic  model.  (In  this  case  a 
nuclear  exchange  model). 

2)  Choose  potential  production  functions. 

3)  Select  experimental  design(s). 

4)  Run  model  for  experimental  design (s). 

5)  Estimate  coefficients  of  production  functions. 

6)  Maximize  the  functions  subject  to  a  cost  con¬ 
straint,  and  analyze  the  results. 

The  rest  of  the  chapter  will  discuss  steps  one  through  five 
in  more  detail.  Chapters  IV  and  V  will  discuss  step  six  for 
the  two  and  three  variable  cases. 

B.  DETERMINISTIC  MODEL 

The  deterministic  model  used  to  generate  data  for  this 
study  is  a  linear  programming  (LP)  model.  This  model  is  a 
strategic  weapons  model.  The  model’s  intent  is  to  maximize 
the  amount  of  damage  against  a  target  base  given  a  set  of 
differing  weapons.  The  LP  model  at  first  glance  is  fairly 
simple.  The  development  of  the  model  though  is  quite  com¬ 
plex  and  includes  weapon  characteristics  (Circular  Error 


Probability,  yield,  reliability,  and  daily  alert  probabil¬ 
ity),  target  characteristics  (hardness  and  vulnerability  to 
particular  kill  mechanism) ,  and  number  of  weapons  assigned 
to  a  target  (one  or  two)  to  yield  a  damage  expected  result 
for  each  weapon,  target,  and  number  of  weapons  used  combina¬ 
tion.  The  model  was  developed  by  Bunnell  and  Takacs  in 
1983.  A  good  description  of  the  model  is  found  in  appendix 
A  of  Graney,  1984. 

A  simplified  example  will  help  to  better  understand 
this  model.  The  example  consists  of  two  weapons  and  two 
types  of  targets.  The  variables  Wft  and  Wg  will  represent 
the  number  of  weapons  of  type  A  and  B  available,  respective¬ 
ly.  The  variables  Ng  and  will  be  the  number  of  targets 
of  type  S  and  type  H,  for  example,  soft  and  hard.  Figure 
3. 1  shows  the  amount  of  damage  expected  (D, ,  D„,  .  .  .  ,  DQ) 
given  a  particular  weapon,  target,  and  number  of  weapons  per 
target  configuration. 

Weapon  A  Weapon  B 

Number  of  Weapons 

per  target  s  1  2  12 

Target  S:  Dj  D 2  Dg 

Target  H*  Dg  D7  Dg 

Figure  3. 1 

The  LP  model  now  consists  of  maximizing  the  total 
destruction  subject  to  limitations  on  the  number  of  weapons 


and  the  number  of  targets.  The  LP  problem  follows  on  the 
next  page.  In  equation  <3.1)f  represents  the  amount  of 

targets  S  attacked  with  weapon  A  using  one  weapon  per 
target,  X2  represents  the  amount  of  targets  S  attacked  with 
weapon  A  using  two  weapons  per  target,  and  so  on.  Equations 
(3.2)  and  (3.3)  are  the  limit  of  targets  S  and  H  available. 
Equations  (3.4)  and  (3.5)  limit  the  amount  of  weapons  avail¬ 
able.  In  these  last  two  equations  the  factor  *2*  appears 
before  the  variables  X2,  X^,  X^,  and  X0  to  account  for  the 

fact  that  a  unit  of  these  variables  is  equivalent  to  two 
weapons  used. 

Maximize 

Z  ■  DiXi+D2X2+D3X3+D4X4+D5X5+D6X6+D7X7+D8XB  (3.1) 

subject  to 
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Admittedly,  the  model  has  some  simplifying  assumptions 
which  do  not  fully  capture  reality.  Theoretically,  a  more 
complex  model  would  have  also  worked  for  the  purpose  of  this 
study.  Unfortunately,  this  was  the  only  unclassified  model 
available.  Still,  the  general  concern  of  this  study  is  to 
examine  production  functions  fitted  to  the  response  surface 


of  a  deterministic  model.  Therefore,  any  LP  model  would 
work.  The  -fact  that  our  model  has  some  potential  military 
application  is  a  bonus  in  this  study. 

The  actual  model  developed  by  Bunnell  and  Takacs  con¬ 
sisted  o-f  ten  target  types  and  -five  weapon  types  (Minuteman 
II,  Minuteman  III,  Peacemaker  (MX),  B-l  bomber,  and  a  sub¬ 
marine  launched  ballistic  missile).  In  order  to  look  -for 
potential  relationships  or  di-f-ferences  among  a  set  o-f  prod¬ 
uction  -functions,  the  model  is  reduced  in  size.  Two  cases 
are  considered,  one  with  two  weapon  types  and  the  second 
with  three  weapon  types.  The  number  of  targets  (not  the 
number  of  target  types)  are  also  reduced,  sixty  percent  in 
the  two  variable  case  and  fifty  percent  in  the  three  vari¬ 
able  case.  The  fifty  and  sixty  percent  reductions  are 
purely  arbitrary.  They  were  meant  to  reduce  the  target 
base,  but  not  so  much  that  the  weapons  could  exhaust  the 
targets  within  the  feasible  range  of  weapons  considered.  In 
other  words,  rather  than  subject  two  or  three  weapons  to  a 
target  base  normally  attacked  with  five  weapons,  a  more 
realistic  assumption  is  to  reduce  the  target  base  and  assume 
the  other  weapons  are  assigned  to  the  rest  of  the  target 
base.  Another  option  is  to  leave  the  other  weapons  in  the 
model  at  some  constant  value,  but  this  causes  difficult 
problems  in  fitting  the  production  functions  to  the  data. 
For  example,  the  weapons  set  at  a  constant  value  do  not 
produce  the  same  effect  for  the  varying  levels  of  the  other 


weapons;  and,  the  LP  model  chooses  different  targets  for  the 
weapons  set  at  constant  level.  If  this  study  had  truly  been 
to  study  the  problem  of  nuclear  exchange  with  the  use  of 
classified  data,  all  weapons  and  targets  would  be  left  in 
the  model,  but  since  its  purpose  is  to  develop  a  methodology 
for  examining  such  decisions  with  the  use  of  other  produc¬ 
tion  functions,  none  of  the  actions  to  the  model  are  of 
significance.  The  data  in  the  model  are  not  classified,  so 
reducing  the  target  base  does  not  affect  the  use  of  the 
model.  If  the  data  were  not  reduced,  the  excess  targets 
cause  the  LP  model  to  place  all  the  weapons  on  one  or  two 
targets,  producing  a  linear  effect  in  the  results.  Reducing 
the  number  of  targets  produces  more  interaction  in  the  model 
between  the  weapons  and  the  targets,  providing  a  problem 
suitable  for  this  research. 

C.  PRODUCTION  FUNCTIONS 

The  functions  or  models  used  in  the  study  to  fit  the 
data  are  classified  into  two  main  groups;  those  character — 
ized  by  constant  elasticity  of  substitution  (CES)  and  those 
with  variable  elasticity  of  substitution  (VES) .  These  func¬ 
tions  have  all  seen  application  in  problems  of  production  in 
economics,  some  much  more  than  others.  The  Cobb-Douglas,  a 
CES  function,  is  seen  widely  in  economic  literature.  On  the 
other  hand,  the  basic  polynomial  (of  order  greater  than  or 
equal  to  two),  a  VES  function,  is  not  seen  very  often  in 
economic  literature  but  appears  commonly  in  other  liters- 


tore,  such  as  operations  research,  as  a  function  for  fitting 
data.  One  purpose  of  this  study  is  to  examine  the  "fitting" 
ability  of  these  various  functions  and  compare  that  Mith  the 
economic  information  gained  from  interpreting  the  coeffi¬ 
cients. 

The  tea  types  of  CES  functions  used  in  this  analysis 
are  the  Cobb-Douglas  and  the  ArroM,  Chenery,  Minhas,  and 
Solow  design.  These  are  referred  to  as  the  Cobb-Douglas 
(C-D)  and  the  CES,  respectively.  The  Arrow,  Chenery, 
Minhas,  and  Solow  function  is  called  the  CES  function  be¬ 
cause  it  is  the  most  general  type  of  CES  function. 

The  C-D  assumes  a  production  process  that  is  somewhat 
"proportional that  is,  the  inputs  are  not  independent  of 
each  other.  One  can  not  have  zero  of  one  input  and  still 
have  output.  The  CES  allows  zero  of  an  input  with  output. 
In  the  LP  model  only  one  weapon  is  needed  for  output; 
therefore,  in  an  aprior  sense  the  C-D  would  not  appear  to  be 
a  goad  model  for  this  problem,  at  least  not  when  one  of  the 
inputs  is  near  zero.  The  CES  function  does  capture  this 
essential  aspect  of  the  problem. 

Four  types  of  VES  functions  are  used  in  the  study,  a 
linear  model  (LIN),  a  polynomial  of  degree  two  with  interac¬ 
tions  (POL),  a  multiplicative  nonhomogeneous  model  (MNH) , 
and  an  additive  nonhomogeneous  model  (ANH) .  The  MNH  produc¬ 
tion  function  was  developed  by  Koenker  and  Perry,  and  the 
ANH  was  developed  by  Sudit,  both  appeared  in  The  Bel 1 


the  Mean  in  econoMic  terns.  The  variance  design  estimated 
polynomial  is  included  because  o-f  its  good  fit  (table  3.1). 

MEAN  ECONOMIC  VALUES 

RATE  OF  TECH. 


MARGINAL 

OUTPUT 

SUBSTITUTION 

PRODUCT 

ELASTICITY 

OF  W„  FOR  W. 

BD  Functions 

W 

"l 

W 

w2 

W 

"l 

W 

2 

L 

*  M  R 

Polynomial 

.82 

.76 

.369 

.566 

0.94 

1.09  1.34 

CES 

.82 

.74 

.371 

.561 

1.10 

1.10  1.10 

Linear 

.82 

.75 

.375 

.574 

1.09 

1.09  1.09 

AklLI 

Hwn 

.85 

.74 

.390 

.565 

0.84 

1.15  1.30 

C-D 

.78 

.75 

.333 

.536 

0.62 

1.04  1.09 

ripen 

.32 

.52 

.  136 

.373 

0.92 

0.61  0.42 

VD  Functions 
Pol ynomi al 

.75 

.69 

.332 

.536 

0.92 

1.09  1.53 

Note:  Mean 

values  are  Weapon  1  = 

225  and 

Weapon 

2  -  375. 

The 

last 

three 

columns  represent  a 

Low  Li 

(50,50), 

Medium,  M:  (225,375),  and  High,  Hi  (400,700),  range 
of  the  inputs. 


TABLE  4.2 


The  first  four  functions  in  table  4.2  agree  quite 
closely  on  some  basic  economic  measures.  As  a  reminder, 
marginal  product  represents  the  marginal  increase  in  output 
for  a  one  unit  increase  in  the  input.  The  output  elastici¬ 
ties  measure  the  percent  change  in  output  for  a  one  percent 
change  in  the  input.  The  rate  of  technical  substitution 
(RTS)  of  *or  is  the  number  of  units  of  which  may  be 
reduced  when  adding  one  unit  of  W^  to  maintain  constant 
output.  In  economic  terms  the  RTS  is  also  the  slope  of  the 
isoquant  for  a  specified  level  of  the  inputs.  An  isoquant 
is  the  combinations  of  inputs  that  maintain  a  particular 
constant  output.  As  can  be  seen  from  the  range  of  the  RTS 
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points  better  than  the  random  data  points.  This  is  not 
unusual  since  the  -functions  were  estimated  from  the  bias 


design  data  points  and  not  the  random  data  points.  The 
polynomial  is  a  remarkedly  good  fit  of  the  points  it  Mas 
estimated  from,  and  still  an  excellent  fit  to  the  random 
points.  Two  anomalities  stand  out  though,  the  additive 
nonhomogeneous,  while  being  a  super  fit  to  the  bias  design, 
is  not  nearly  as  good  a  fit  to  the  random  points.  On  the 
other  hand,  the  CES  is  the  third  best  fitting  fuction  to  the 
bias  design  and  the  best  fitting  function  to  the  random  set. 

In  terms  of  the  two  measures  of  fitting  ability  (SDS 
and  percent  fit),  the  functions  were  subjectively  broken 
into  three  major  groups  which  could  be  called  “marginal , " 
“better,"  and  “best."  The  “marginal"  group  includes  the  MNH 
and  C-D.  The  “better"  group  contains  the  Linear  and  the 
ANH.  The  "best"  group  is  the  CES  and  the  polynomial. 

B.  INTERPRETATION  VALUE 

Another  method  for  determining  which  production  func¬ 
tion  to  use  is  to  look  at  its  economic  interpretation.  In 
the  following  sections  of  this  chapter,  each  of  these 
functions  will  be  analyzed  with  respect  to  its  economic 
interpretation  and  its  subsequent  implications.  Before  pro¬ 
ceeding  into  the  individual  functions,  it  would  be  of  inter — 
est  to  see  how  the  functions  compare  to  each  other  in  eco¬ 
nomic  terms.  Table  4.2  compares  the  production  functions  at 
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functions  fit  extremely  well.  Unfortunately,  it  is  not  so 
clear  which  function  fits  best  because  the  actual  response 
is  Measured  in  neither  of  these  extremes.  The  LP  output 
ranged  from  zero  to  one  thousand  with  the  average  response 
around  five  hundred.  The  SDS  values  in  table  4.1  can  be 
used  to  compare  the  functions  against  each  other  by  exam¬ 
ining  the  difference  in  the  magnitudes  of  error  between  the 
function. 

Another  possible  measure  of  fit  is  the  percent  fit  of 
the  predicted  responses  to  the  actual  responses.  Here,  one 
must  decide  what  is  an  acceptable  percentage.  If  that  value 
is  defined  as  95  percent,  then  three  of  the  functions  meet 
that  criteria  as  well  as  one  of  the  other  functions  esti¬ 
mated  using  the  variance  design. 

Finally,  it  ought  to  be  mentioned  that  the  different 
measures  of  fit  could  give  conflicting  results.  In  which 
case  the  acceptance  of  a  function  may  more  depend  on  its 
interpretational  value  than  on  whether  it  is  a  slighty 
better  fit  than  another  function.  It  is  interesting  to  note 
in  table  4.1  that  the  percent  fit  of  the  polynomial  is  less 
than  that  of  the  ANH,  yet  the  SDS’s  suggest  that  the  poly¬ 
nomial  would  fit  better. 

Examining  table  4.1,  one  sees  that  each  of  the  six 
functions  give  slightly  different  conclusions  when  applied 
to  the  bias  data  base  as  opposed  to  the  random  data  base. 
As  one  would  expect,  the  functions  fit  the  bias  design  data 
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fitting  function.  Sines  all  the  functions  are  estimated 
from  the  same  design,  the  Sum  of  the  Deviations  Squared 
(SDS)  of  the  bias  design  need  only  be  compared.  Table  4.1 
compares  the  actual  SDS  of  the  functions  estimated  with  the 
bias  design  when  applied  to  the  all  bias  design  and  the  SDS 
when  applied  to  the  random  data  base.  The  last  column  in 
table  4.1  is  the  average  for  all  the  points  of  one  minus  the 
percent  error  (the  predicted  response  minus  the  actual  re¬ 
sponse,  divided  by  the  actual  response) . 


TWO  Variable  Functions  Estimated  Using  Bias 

Desi gn 

****DATA 

POINTS**** 

PERCENT 

BIAS 

RANDOM 

FIT  TO 

SDS 

SDS 

RANDOM 

Linear 

453. 

1243. 

.94 

Polynomial 

50. 

555. 

.96 

Cobb -Doug 1 as 

4900. 

26974. 

.91 

Mult.  Non Horn. 

1242. 

10305. 

.88 

Add.  NonHom. 

89. 

1392. 

.98 

CES 

349. 

487. 

.98 

Notes  BIAS  design 

contains 

nine  points. 

Random  design  contains  fifteen  points. 
TABLE  4.1 


Before  examining  table  4.1,  a  couple  of  points  need  to 
be  made.  Alone,  the  SDS  does  not  tell  how  good  a  fit  the 
functions  are  to  the  surface.  For  instance,  if  the  SDS  is 
the  only  measure  of  fit  and  if  the  response  is  measured  in 
units  no  greater  than  10,  the  results  of  table  4.1  would 
suggest  that  none  of  the  functions  fit  well.  On  the  other 
hand,  if  the  response  is  measured  in  billions,  all  the 
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A-  PREDICTION  VALUE; 

There  arc  two  important  elements  in  choosing  a  function 
to  best  represent  the  response  of  a  surface.  First,  it 
should  be  able  to  accurately  predict  a  response  for  a  given 
set  of  inputs,  and  secondly,  it  ought  to  make  sense  in  the 
interpretation  of  the  coefficients,  that  is,  either  confirm 
an  aprior  theory  or  give  insight  that  is  not  contradictory 
to  the  process  generating  the  response.  In  the  first  part 
of  this  chapter,  prediction  will  be  the  measure  of  effec¬ 
tiveness  used  in  evaluating  the  production  functions. 

As  was  discussed  in  the  previous  chapter,  in  order  to 
compare  the  various  functions,  a  random  set  of  inputs  to  the 
real  surface  was  generated  and  input  for  each  of  the  produc¬ 
tion  functions.  The  random  data  set  was  not  used  to  esti¬ 
mate  any  functions,  but  rather  to  compare  the  predicted 
values  of  the  estimated  functions  against  each  other.  The 
benefit  of  the  random  data  set  is  that  all  the  functions  are 
compared  against  a  common  data  set.  Since  some  of  the 
functions  were  estimated  by  different  data  sets  (variance 

design  or  bias  design),  there  was  no  common  basis  to  compare 
the  functions.  The  results  (Table  3.1)  showed  that  in  all 
cases  the  bias  design  provided  as  good  a  fit  or  better  to 
the  surface  than  the  variance  design  employed.  Therefore, 
the  functions  with  coefficients  estimated  with  the  bias 
design  are  the  only  functions  to  be  considered  for  the  best 
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E.  eSBItfi  CQ§I  IQ  IHg  PROBLEM 


At  the  start  of  the  study,  some  approaches  were  dis¬ 
cussed  for  adding  cost  into  the  analysis.  One  suggestion 
was  adding  a  cost  variable  in  the  -function  estimated  using  a 
response  surface.  In  effect  this  would  have  meant  adding 
another  constraint,  a  cost  constraint,  to  the  LP  model. 
Unfortunately,  this  does  not  work  well  since  it  produces 
quite  conflicting  results.  The  output,  damage  produced, 
really  depends  on  the  number  of  weapons  or  the  total  dol¬ 
lars.  It  doesn’t  depend  on  any  combination  of  these  two 
factors.  Either  there  is  sufficient  funding  and  weapons 
restrain  the  production  or  there  are  sufficient  weapons  and 
funding  restrains  the  problem.  Inputting  a  funding  level 
and  weapon  level  into  an  equation  does  not  tell  which  is  the 
real  driver  in  production  (damage),  and  the  response  surface 
would  tend  to  average  the  effect  of  the  competing  con¬ 
straints. 

A  more  realistic  approach  to  adding  cost  to  the  problem 
is  to  maximize  the  particular  production  function  subject  to 
a  cost  constraint.  This  was  done  in  the  two  variable  case 
using  a  Lagrangian  technique.  For  more  than  two  variables, 
a  nonlinear  computer  optimization  program  is  really  needed. 
The  results  of  adding  cost  have  quite  different  implications 
for  the  varying  production  functions  and  will  be  discussed 
more  in  the  next  chapter. 


sign  the  all  bias  design  would  require  seventy-six  design 
points,  the  variance  design,  only  f ifty— Four . 

In  dealing  with  the  problem  of  multicol  1  inearity  in  the 
all  bias  designs  a  couple  of  simple  rules  are  used  in  the 
selection  of  variables  to  include  in  the  model.  If  during 
the  stepwise  regression  the  adjusted  R  square  decreases  or 
the  coefficients  change  drastically  in  sign,  then  the  vari¬ 
able  is  not  included  and  the  process  is  stopped  there.  In 
cases  where  a  coefficient  was  not  estimated,  it  was  usually 
the  last  variable  left  to  be  included  in  the  regression. 

The  Box  and  Bhenken  design  used  in  this  study  requires 
the  use  of  fifteen  design  points  for  the  three  variable 
case.  These  fifteen  points  include  a  center  point  repli¬ 
cated  three  times.  For  a  deterministic  model  there  is  no 
need  to  replicate  the  center  point  since  one  will  get  the 
same  response  each  time.  An  alternate  idea  to  leaving  the 
replications  out  is  to  use  psuedo-repl i cations,  that  is, 
moving  off  the  center  one  unit  in  any  direction.  An  advan¬ 
tage  to  doing  this  is  that  it  provides  more  "weight"  to  the 
center  of  the  design.  Also,  since  the  all  bias  design 
requires  fifteen  design  points,  it  makes  the  two  designs  at 
least  equal  in  number  of  design  points.  There  is  no  notic- 
able  change  in  the  orthogonality  of  the  variance  design  by 
the  use  of  psuedo  replications. 


The  first  drawback  is  the  lack  of  orthogonality  in  the 
design  for  models  of  order  two.  The  variance  design  pos¬ 
sesses  a  high  degree  of  orthogonality.  Because  of  this  fact 
and  the  lack  of  variance  error  in  the  data,  all  variables 
can  be  included  when  estimating  using  multiple  regression. 
This  is  not  always  the  case  with  the  bias  design.  Because 
of  the  lack  of  orthogonal i ty ,  there  is  a  problem  with  multi- 
collinearity.  The  coefficients  can  change  drastically  if 
too  many  variables  are  added  to  the  model.  A  possible 
solution  would  be  to  use  an  orthogonal  all  bias  design  based 
somewhat  on  the  Box-Bhenken  designs,  or  on  some  other  or — 
thogonality  technique. 

Another  benefit  of  orthogonal  designs  is  that  the  error 
contribution  of  each  term  can  be  easily  determined.  This 
allows  one  to  see  which  terms  best  fit  the  surface  and  which 
terms  are  most  important  in  predicting.  This  can  be 
accomplished  because  the  variance/covariance  matrix  of  an 
orthogonal  design  is  zero  everywhere  but  along  the  diagonal. 
The  diagonal  then  represents  the  contribution  of  error  from 
each  term  in  the  design. 

The  second  drawback  with  the  all  bias  design  is  the 
number  of  points  to  estimate.  In  the  two  cases  in  this 
study  the  lumber  of  designs  points  are  the  same,  but  in  any 
design  for  more  than  four  variables  the  all  bias  design 
would  require  the  calculation  of  more  data  points  than  the 
Box-Bhenken  designs.  For  example,  for  a  five  variable  de- 


TWO  Variable  Function  Estimated  By 
Variance  Design  Bias  Design 
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SDS 

PRODUCTION  FUNCTION 

Linear 

87818. 

(.897) 

86875. 

(.899) 

Polynomial 

6014. 

(.981) 

4187. 

(.980) 

Cobb-Douglas 

1095204. 

(.608) 

228665. 

(.876) 

Mult.  Non Horn. 

777413. 

(.671) 

128150. 

(.909) 

Add.  Non Horn. 

+ 

+ 

CES 

+ 

45706. 

(.932) 

—  not  estimated 
Table  3.1 


The  analysis  of  the  results  of  table  3.1  as  they  apply 
to  the  best  function  to  use  and  its  interpretation  will  be 
more  fully  discussed  in  Chapters  IV  and  V.  What  is  inter¬ 
esting  to  note  here  is  the  improvement  in  accuracy  using  a 
bias  design  against  a  more  conventional  variance  design.  It 
should  not  be  forgotten  that  the  data  have  no  variance  since 
it  came  from  a  deterministic  model  which  implies  that  all 
the  error  is  bias  error.  The  MNH  and  ANH  functions  bene- 
fitted  the  most  with  the  use  of  the  bias  design,  at  least  in 
the  two  variable  case.  Unf ortunately,  there  are  some  draw¬ 


backs  in  the  use  of  an  all  bias  designs  to  deterministic 


Multiple  regression  is  used  to  -fit  the  experimental 


design  responses  to  the  production  functions.  The  CES,  as 
mentioned  above,  required  the  use  of  a  nonlinear  estimation 
technique. 

Although  it  is  not  the  specific  intention  of  this 
study,  a  subgoal  was  to  determine  the  practicality  of  using 
all  bias  designs  to  fit  a  function  to  data  of  a  determini¬ 
stic  model.  To  measure  the  various  fitting  abilities  of 
the  functions,  those  estimated  using  a  VD  and  those  using  a 
BD,  a  random  set  of  data  points  to  the  real  surface  is 
generated  and  applied  to  each  estimated  function.  The 
error  is  computed  as  the  expected  response  minus  the  actual 
response.  The  Sum  of  the  Deviations  Squared  (SDS>  is  com¬ 
pared  for  each  of  the  functions  estimated  with  different 
designs  (VD  or  BD) .  In  the  application  of  variance  designs 
and  bias  designs  to  fitting  these  functions  it  was  dis¬ 
covered  that  the  bias  designs  provide  a  better  fit.  Table 
3.1  shows  the  result  of  this  excursion  using  15  randomly 
generated  data  points  for  the  two  variable  case  and  the 
results  of  25  random  points  for  the  three  variable  case. 
The  number  in  the  parenthesis  could  be  called  the  "percent 
fit  to  the  surface."  It  is  one  minus  the  average  of  the 
percent  error.  The  percent  error  is  the  absolute  value  of 


the  predicted  response  minus  the  actual  response  divided  by 
the  actual  response. 


All  Bias  Designs 
Variablesi  1  2  ...  k 
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Note:  g  “  (CiiD)  and  a  *  2  g  -for  rotatability 

where  k  >  #  of  variables  and  Ciil  *  l/(k+2>. 

Figure  3.3 

and  the  Maximum  value  is  set  equal  to  positive  one.  The 
mean  or  middle  value  is  then  zero.  The  coded  values  can  be 
used  directly  in  the  regression  analysis  or  the  actual 
variable  values.  The  range  used  in  this  analysis  is  0-450 
for  weapon  one,  0-730  for  weapon  two,  and  0-1040  for  weapon 
three. 

Some  of  the  functions  estimated  are  not  second  order. 
The  C-D  and  linear  function  are  first  order  functions.  In 
the  case  of  the  ANH,  it  is  not  clear  what  type  of  design  the 
function  would  require.  The  second  order  design  is  the  best 
choice  for  a  design  since  a  second  order  design  will  esti¬ 
mate  first  order  and  second  order  models  correctly.  The  CES 
model  cannot  be  estimated  using  multiple  regression  and  a 
nonlinear  approach  using  the  same  experimental  design  was 
employed  with  the  CES. 


The  -first  set  o-f  designs,  the  variance  designs  <VD> , 
are  a  -full  -factorial  for  the  two  variable  case  and  a  rotat¬ 


able  second  order  design  developed  by  Box  and  Behnken  for 
the  three  variable  case  (Box  and  Bhenken,  1960x455).  Both 
designs  possess  a  high  degree  of  orthogonality.  The  designs 
for  the  two  and  three  variable  case  are  given  in  figure  3.2. 
The  Box  and  Bhenken  work  contains  designs  for  problems  of  up 
to  sixteen  variables. 


DESIGN  POINTS 
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Figure  3.2 


The  second  set  of  designs,  all  bias  designs  (BD) ,  seek 
primarily  to  minimize  bias.  These  designs  are  developed 
from  the  work  of  Box  and  Draper  summarized  in  Myers*  book  on 
RSM  (Box  and  Draper,  1959:622-654).  They  require  the  cal¬ 
culation  of  ’a’  and  ’g’  from  the  second  moment,  generally 
noted  in  the  literature  as  *Cii3’.  The  general  design  is 
given  in  figure  3.3. 

In  both  figure  3.2  and  3.3,  there  is  an  assumption  that 
coding  is  being  used  in  the  design.  For  these  designs  it  is 
assumed  that  the  mininum  value  is  set  equal  to  negative  one 


times  751)  for  the  two  variable  case  and  3. 5258774E+08  runs 
(451  times  751  times  1041)  for  the  three  variable  case.  In 
the  three  variable  case,  with  each  run  taking  0.40  CPU 
seconds,  this  would  amount  to  four  and  one-half  years  of 
computer  run  time.  Obviously  the  resource  cast  is  too  great 
to  completely  estimate  the  function  by  data  points,  and  data 
points  alone  do  not  tell  much  about  the  relationships  be¬ 
tween  the  weapons.  For  this  reason  specific  points  are 
chosen  and  measured  and  functions  are  fitted  to  these  points 
to  estimate  the  true  surface.  Experimental  designs  provide 
the  vehicle  to  estimate  the  surface  with  a  mininum  of  data 
points. 

Most  designs  used  in  fitting  a  function  to  the  response 
attempt  to  minimize  variance.  In  fitting  a  function  to  the 
response  of  a  deterministic  model  the  only  error  involved  is 
bias  error;  there  is  no  variance  error.  Bias  measures  the 
error  in  choosing  the  wrong  function  to  fit  the  surface.  Ir* 
choosing  an  experimental  design  it  seemed  appropriate  to  use 
a  design  which  minimized  bias.  Actually,  Myers  concluded 
from  the  work  of  Box  and  Draper  that  designs  which  minimize 
variance  do  a  very  good  job  at  minimizing  bias  also  (Myers, 
1976*208,  230).  Still,  in  this  study  two  experimental  de¬ 
signs  are  used,  a  variance  minimizing  design  and  a  bias 
minimizing  design.  Both  are  second  order  designs,  that  is, 
designs  specifically  for  fitting  a  second  order  equation  to 


the  data 


Journal  of  Economics  and  Management  Sc; lance.  The  linear 
model  is  the  simplest  of  the  four,  while  the  polynomial  is 
the  most  complex.  The  polynomial  model  is  used  extensively 
in  non-economic  research.  There  is  good  reason  for  this 
wide  use  of  the  polynomial,  because  it  is  a  quasi -Taylor 
expansion  of  the  true  and  unknown  function  (Box  and  Draper, 
1971:734). 

Of  these  four  functions  the  IWH  would  least  likely  fit 
the  production  process  represented  in  the  LP  model.  The  MNH 
requires  positive  values  of  all  inputs  before  there  is  any 
output.  Still  the  MNH  is  chosen  for  use  to  see  how  it  fits 
away  from  the  origin  and  to  see  what  kind  of  economic  in¬ 
sight  it  gives. 

All  of  these  models  have  been  discussed  in  more  detail 
in  chapter  two  of  this  report  and  further  questions  should 
be  directed  there  or  to  the  references  in  the  bibliography. 

D.  EXPERIMENTAL  DESIGN 

The  selection  of  an  experimental  design  proved  to  be 
one  of  the  more  difficult  aspects  of  this  study.  Since  the 
data  do  not  come  from  an  actual  experiment,  but  from  a 
deterministic  model,  there  is  no  variance  error  in  the  data. 
The  data  represent  the  actual  or  true  surface.  In  fact  if 
one  had  the  computer  resources  he  could  calculate  the  re¬ 
sponse  for  every  input  combination.  In  the  simple  problem 
used  in  this  study  that  would  amount  to  338,701  runs  (4S1 


values,  the  slope  of  the  isoquants  do  not  change  much  for 
differing  levels  of  the  inputs.  These  results  also  seen  to 
suggest  that  the  MNH  is  not  a  good  model  of  the  response 
surface.  It’s  economic  interpretation  is  completely  out  of 
line  with  the  other  functions  which  also  fit  the  surface 
better.  The  Cobb-Douglas,  while  not  agreeing  with  the  other 
functions  as  well,  is  not  that  poor.  Additionally,  there  is 
evidence  to  suggest  the  C-D*s  fit  is  fairly  good  away  from 
the  origin.  This  will  be  discussed  more  fully  in  the  sec¬ 
tion  on  the  Cobb-Douglas. 

In  choosing  a  best  fit,  these  last  two  sections  have 
examined  both  the  predicting  ability  and  a  general  economic 
interpretation  at  the  mean.  In  order  to  more  closely  ex¬ 
amine  the  economic  interpretation,  each  function  needs  to  be 
looked  at  in  greater  detail.  The  following  sections  will 
examine  the  functions  beginning  with  the  better  fitting  and 
ending  with  the  least  best  fit. 

1.  IHE  CES  PRODUCTION  FUNCTION.  The  CES  production 
function  is  a  very  common  model  in  economic  literature. 
Each  of  the  parameters  estimated  in  the  CES  function  have 
very  descript  properties.  The  form  of  the  CES  is  as 
follows: 

-p  -p  -v/p 

Y  -  A<«K  +  <1— «>L  > 

A  is  an  efficiency  parameter  (A>0);  P  is  a  substitution 
parameter  <-l<P<inf inity) ;  «  is  a  distribution  parameter 
(v<a£l);  V  measures  returns  to  scale  <V>0) . 


The  CES  -function  was  estimated  with  the  bias  design 


using  a  nonlinear  approach.  Initially  none  o-f  the  values 
were  restricted.  Below  is  the  result  when  the  parameters 
are  allowed  to  range  -freely. 

Y  *  2. 23 (  . 53W1 *  07  +  .47W1*07)* 95/1,07 

Surprisingly,  one  o-f  the  parameters  when  estimated  did 
not  fall  in  the  prescribed  range.  The  P  parameter  which 
measures  the  ability  to  substitute  among  the  inputs  was  just 
barely  less  than  negative  one  (-1.07).  A  P  value  of  nega¬ 
tive  one  indicates  perfect  ease  of  substitution.  A  large  P 
value  would  indicate  a  proportional  process.  For  example, 
mowing  lawns  is  a  proportional  process  with  substitution 
extremely  difficult.  If  there  are  three  lawnmowers  and 
three  individuals  to  mow  the  lawns,  one  cannot  substitute 
lawnmowers  for  more  individuals  to  mow  the  lawns  without 
losing  some  production.  Perfect  ease  of  substitution  im¬ 
plies  the  opposite;  it  is  very  easy  to  substitute  among  the 
inputs  and  maintain  production. 

Perfect  ease  of  substition  is  a  seemingly  obvious  re¬ 
sult  for  the  case  of  two  weapons  producing  damage.  There¬ 
fore,  P  was  set  equal  to  negative  one  and  estimated  again  by 
nonlinear  estimation  techniques. 


Y  *  2.34 (  .52Wt  +  ,4BW2> 


(4.1) 


One  interpretation  of  equation  4.1  is  that  the  process 
is  slightly  more  weapon  one  intensive  than  weapon  two,  but 
what  does  this  mean?  In  a  more  typical  capital-labor  pro- 
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duction  problem,  if  capital  had  a  higher  intensity  value 
than  labor,  it  means  that  the  technology  is  more  capital 
intensive.  It  appears  that  in  this  context,  the  intensities 
measure  the  capabilities  of  the  weapons.  The  ratio  of  the 
intensities  is  the  slope  of  the  isoquant,  for  all  isoquants. 
This  will  became  clearer  later  in  the  discussion,  hence,  the 
ratio  of  the  intensities  is  the  rate  of  technical  substitu¬ 
tion.  Therefore,  the  intensities  reflect  the  mean  capabili¬ 
ties  of  the  weapons. 

Also,  equation  4.1  shows  that  the  process  has  dimin¬ 
ishing  returns  to  scale  (.94  <  1).  There  is  proportionately 
less  payoff  for  more  weapons.  This  makes  sense  since  as  one 
increases  the  number  of  weapons,  he  must  attack  progres¬ 
sively  more  difficult  targets.  After  the  easier  targets 
have  been  attacked  there  are  only  incremental  increases  in 
destruction  for  more  weapons. 


The  efficiency  parameter  in  equation  4. 1  is  a  more 
difficult  parameter  to  interpret  unless  there  is  another 


production  process  to  compare  it  with. 

Figure  4.1  is  a  plot  of  equation  4.1.  The  X-axis 
(Weapon  1)  ranges  from  0  to  450}  the  Y-axis  (Weapon  2) 
ranges  from  0  to  750.  The  plot  looks  very  much  like  a 
plane,  yet  fits  much  better  than  the  linear  production 
function  which  is  a  plane. 

The  most  interesting  aspect  of  these  varying  production 
functions  is  the  different  interpretations  one  gets  when 
cost  is  included  in  the  analysis.  To  identify  this  rela¬ 
tionship,  the  production  function  will  be  maximized  subject 
to  a  cost  constraint. 

-P  -p  — v/P 

Maximixe  Y  «*  A(«W^  +  (1-«)W2 > 

Subject  to  TC  -  ♦  C2W2 

Notei  TC  is  total  funding,  is  the  price  of  weapon  one  and 
C2  is  the  price  of  weapon  two. 

OR 

Max  L  -  A(«W~P  +  (1-«)W2P)"V/P  -  +  C2W2  -  TC)  (4.2) 

Equation  4.2  is  the  Lagrangian.  Next,  the  partials  are 
taken  and  set  equal  to  zero. 


(— V/P) A(«W~P  +  (1-«)W^P)  <’V/P)~1(-P)«W7P”1  -  jiCj  -  0 

(4.3a) 


(-V/P)A(«W^P 


(1-«)W^P)  <”V/P)_1  (-P)  (l-«)W^P"'1-pC2  -  0 

(4.3b) 


—  -  C1W1  +  C2W2  -  TC 


(4.3c) 


Solving  -for  >4  in  the  first  two  equations  and  then  setting 
them  equal  to  each  other  yields  the  -following  relationship: 


=  c«c2/(i-«)c1i1/<p+1,w2 


(4.4) 


Substituting  equation  4.4  into  equation  4.3c  gives  the 
amount  of  W2  required  to  maximize  damage,  the  demand  equa¬ 
tion  (equation  4.5). 


c1c«c2/(i-«)c1i1/(p+l>  +  c2 


(4.5) 


The  substitution  coefficient  is  c  »  1/(P+1).  Since 
P  *  -1,  there  is  absolute  ease  of  substitution  (<r  *  infini¬ 
ty).  Within  the  demand  equation  4.5  is  a  very  important 
ratio  (equation  4.6).  Since  e  is  infinite,  there  are  three 
potential  results.  If  the  value  within  the  brackets  of 
equation  4.6  is  equal  to  one,  then  equation  4.5  simplifies 
to  the  form  W2  =  TC/CC^  +  C^,  and  one  buys  an  equal  number 
of  weapons  one  and  two.  If  the  ratio  within  the  brackets  is 
less  than  one,  then  equation  4.5  simplifies  to  W2  =  TC/C2 
and  only  weapon  two  is  bought.  If  the  ratio  within  the 
brackets  is  greater  than  one,  the  denominator  goes  to  infin¬ 
ity  causing  no  weapon  two’s  to  be  bought,  only  weapon  one. 


C«C2/(1-«)C11 


(4.6) 


The  only  case  where  the  ratio  can  be  equal  to  one  is  if  the 
following  holds: 


«/(!-«)  »  C./C_ 


(4.7) 


Substituting  the  actual  values  implies  that  in  order  to 
buy  a  mix  of  weapons  the  ratio  of  the  cost  of  weapon  one  to 
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other , 
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result  is  gained.  It  is  clear  then,  that  the  ratio  of  the 
intensities  is  the  rate  of  technical  substitution. 

Since  the  ratio  of  the  intensity  parameters,  a/l-a,  is 
the  slope  of  the  isoquant,  and  the  ratio  of  prices,  C^/C^t 
is  the  slope  of  the  budget  curve,  for  the  case  where  the 
ratio  of  the  prices  is  equal  to  the  ratio  of  the  intensi¬ 
ties,  any  combination  of  weapon  one  and  weapon  two  is  the 
best  economic  choice.  The  math  suggested  an  equal  amount  of 
each  weapon,  but  for  this  special  case  there  are  an  infinite 
number  of  solutions.  The  decision  to  buy  all  of  either 
weapon  is  one  of  these  infinite  number  of  solutions.  Hence, 
for  any  ratio  of  prices,  the  decision  to  buy  all  of  one 
weapon  is  the  best  economic  choice. 

The  conclusion,  that  in  all  cases,  one  only  buys  one 
type  of  weapon  and  not  a  mix,  does  not  seem  to  be  like  one's 
intuition  of  a  strategic  weapons  problem.  Unless  both  weap¬ 
ons  possess  similar  capabilities  across  a  range  of  different 
targets,  one  would  not  expect  this  answer.  The  weapons  are 
generic,  and  no  attempt  was  made  to  identify  what  a  weapon 
most  closely  resembles  since  the  data  are  only  unclassified 
estimates.  This  conclusion  hinges  on  the  fact  that  the 


weapons  are  perfect  substitutes  in  economic  terms 


The 


application  is  that  cost  effective  weapon  decisions  can  be 
made  based  only  on  the  price  of  the  weapons  and  the 
intensity  parameters  of  the  CES  production  function. 


2.  JHE  POLYNOMIAL  PRODUCTION  FUNCTION.  The  polynomial 
is  probably  the  most  obscure  of  the  production  functions 
estimated  in  this  study.  The  coefficients  of  the  polynomial 
do  not  have  the  descript  properties  of  the  CES,  and  it  takes 
some  work  to  find  its  economic  interpretation.  The  form  of 
the  polynomial  for  two  variables  is  given  in  equation  4.8. 

V  -  eO  +  'lWl  +  *2*2  +  *12W12  ♦  '11“?  +  *22W2  <4-8> 

The  0O  term  implies  some  ability  to  produce  when  all 

the  inputs  are  zero.  Obviously,  this  does  not  make  economic 

sense.  The  function  can  be  estimated  without  the  pQ  term, 

but  the  fit  is  a  little  better  with  it.  Since  it  is  not 

very  large,  the  pQ  term  was  left  in  the  model.  Equation  4.9 

and  4.10  are  the  results  of  estimating  the  polynomial 

through  multiple  regression  for  the  two  variable  case. 

Figure  4.2  is  a  three  dimensional  plot  of  equation  4.9. 

Bias  Design  Estimated 

Y  -  -21.8  +  .91*^  +  .99W2  -  2.  4  <  lO*"4)  Wj  k»2  -  2.4(10“4>W^ 

<4.  9) 


Variance  Design  Estimated 

Y  -  -12.4  +  .96Wl  +1.0AW2  -  3.9(10“4) -  1.4(10”4)W^ 

-  3.8<10"4)W^  <4. 10) 
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UEAPON  TWO  < 0-750) 

FI CUBE  4.2 

The  polynomial  is  not  homogeneous  and  therefore,  the 
returns  to  scale  cannot  be  measured.  The  elasticities  of 
output  can  be  determined  though,  using  the  standard  defini¬ 
tions 

Ey(X  -  <X/Y> <bY/bX) 

It  appears  that  because  of  scale  and  the  non— homogeneous 
nature  of  the  function,  the  output  elasticities  change  for 
differing  levels  of  production.  Table  4.3  is  the  elastici¬ 
ties  at  the  different  extreme  combinations  of  the  inputs. 
The  elasticities  agree  closely  for  both  designs.  The  bias 


design,  in  general,  has  slightly  larger  values. 


INPUT 

BIAS 

DESIGN 

VARIANCE 

DESIGN 

LEVEL 

ELAST 

ELAST 

ELAST 

ELAST 

W1 

W2 

Ml 

M2 

SUM 

Ml 

M2 

SUM 

SO 

50 

.62 

.66 

1.28 

.54 

.58  1 

.12 

50 

700 

.06 

.76 

.82 

.06 

.62 

.68 

400 

50 

.93 

.11 

1.04 

.84 

.11 

.95 

400 

700 

.35 

.46 

.91 

.29 

.33 

.62 

MEANS  225 

375 

.37 

.57 

.94 

.33 

.51 

.84 

TABLE  4.3 


If  the  polynomial  was  homogeneous,  the  elasticities 


would  not  be  affected  by  scale  changes.  The  Cobb-Douglas, 
■for  example,  has  elasticities  o-f  output  that  do  not  change 
for  different  levels  of  output.  The  CES  elasticities  of 
output  do  change  at  the  varying  levels  of  production,  but 
always  sum  to  the  same  value.  As  can  be  seen  from  table 
4.3,  the  polynomial's  elasticities  do  change  and  do  not  sum 
to  the  same  value. 

The  ease  of  substitution  is  a  difficult  parameter  to 
exact.  Equation  4.11  is  the  formula  used  in  estimating  the 
elasticity  of  substitution  (McFadden,  1963:75).  One  would 
expect  it  to  be  high  given  a  knowledge  of  the  production 
process  (weapons  producing  damage),  yet  it  is  not  infinite 
(implying  perfect  ease  of  substitution)  as  in  the  CES  case. 


«.  »  (1/x.f.  +  1/x  f  )/(-f../f?  +  2f.  /f.f  .  -f  /f2> 

i  j  i  i  j  i  ui  i j  i  j  Jj  3 


(i#j) 


(4.11) 


The  mean  value  for  i*  one  hundred  twenty  three;  it 
ranged  from  fifteen  to  fifteen  hundred.  The  low  values 
occur ed  at  the  maximum  levels  of  the  weapons. 

Figure  4.3  is  a  plot  of  four  isoquants.  An  isoquant  is 
the  combination  of  inputs  that  maintain  a  constant  produc¬ 
tion.  As  can  be  seen  from  the  figure,  the  isoquants  are 
fairly  flat,  especially  within  the  actual  range  of  the 
weapons.  An  isoquant  that  is  a  straight  line  implies  per¬ 
fect  ease  of  substitution.  This  figure  supports  the  CES 
result  of  perfect  ease  of  substitution. 


FIGURE  4.3 


Maximizing  the  polynomial  subject  to  cost  gives  the 
Lagrangian  (equation  4.12). 

Max  L  =  +  p1WJ  +  P2W2  +  ^12W1W2  +  ■*nWi  +  <®22W2 

-  XtCjWj  +  C2W2  -  TC)  (4.12) 
Taking  the  partial  derivatives  yields  the  following  set  of 
equations: 

bL 

—  *  #»x  ♦  #*i2W2  +  2pllWl  “  wCl  “  0  (4.13a) 

Wt 

bL 

~  ”  P2  +  ^12Wt  +  2^22W2  "  Wc2  *  0  <4. 13b) 

W2 
bL 

—  *  C  W  +  C2W  -  TC  =  0  (4.13c) 

bp 

Eliminating  q  in  equations  4.13a  and  4.13b  gives  the 
following  relationship: 


Substituting  into  4.13c  gives  the  relationship  that  tells 
how  aany  of  either  weapon  one  buys  depending  on  the  total 


funding,  TC,  and  prices,  . 


TC(2*11C2  -  P12V  -  p2q  ♦  »1C2C1 

2<#,22C1  +  PilC2  ~  012C2C1> 


(4.15) 


It  is  not  clear  from  looking  at  equation  4.15  what  will 
happen  when  either  total  cost  or  an  individual  weapon  cost 
changes.  As  an  exercise,  these  equations,  with  the  values 
from  the  polynomial,  were  tested  on  some  fictional  data. 
The  use  of  fictional  cost  creates  no  problem  since  it  is 
only  measuring  the  relative  effect  of  costs.  Table  4.4 
illustrates  the  case  where  the  cost  of  the  weapons  are  held 
the  same,  but  there  is  an  increase  in  funding. 


FUNDINGS  4000 

Cl* 

5.0 

C_s 

4.0 

BD-POLYNOMIAL 

wjs 

244 

W?s 

695 

VD-POLYNOMIAL 

u*. 

wr 

388 

Ws 

2 

515 

FUNDINGS  4200 

ci* 

5.0 

C,s 

4.0 

BD-POLYNOMIAL 

wjs 

362 

W?s 

597 

VD-POLYNOMIAL 

445 

w2. 

W  . 

494 

TABLE  4.4 

It  is  interesting  to  note  that  the  number  of  weapons  of 
type  two  decrease  as  funding  is  increased.  In  economic 
terms,  this  suggest  an  inferior  good.  An  example  of  an 
inferior  good  might  be  potatoes  in  Ireland  in  the  eighteenth 
century.  Potatoes  were  needed  for  survival,  they  were  a 
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PERCENT  OF  WEAPON  LEU EL 

FIGURE  4.4 


basic  necessity.  Yet,  if  income  went  up,  people  would  buy 
less  potatoes  for  food  and  more  of  other  staples. 

Is  weapon  two  really  an  inferior  good  when  compared 
with  weapon  one?  A  clear  measure  of  the  capability  of  a 
weapon  is  its  marginal  product.  The  marginal  products  are 
easily  computed;  they  are  simply  the  first  partials  of  the 
function.  Figure  4.4  is  a  graph  of  the  first  partials  over 
the  range  of  the  weapons.  Percent  of  weapon  level  refers  to 
the  percentage  of  weapons  used  within  the  limits  of  the 
weapon  range.  For  example,  100  percent  usage  of  weapon  one 
is  450  weapons;  100  percent  of  weapon  two  is  750  weapons. 

As  the  level  of  production  increases,  weapon  one 
becomes  a  better  weapon  than  weapon  two.  What  is  happening 
in  the  model  is  that  weapon  two  represents  a  somewhat  in¬ 
accurate  but  large  megatonnage  weapon.  For  soft  targets, 
accuracy  is  not  very  important,  and  a  larger  payload  (mega- 
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tonnage)  achieves  more  destruction.  As  the  level  of  produc¬ 
tion  increases  and  more  difficult,  harder  targets  are 
attacked,  weapon  one,  which  represents  good  accuracy  and 
seall  payload,  performs  better.  Hard  targets  are  more  sus¬ 
ceptible  to  accurate  weapons  than  large  payload  weapons. 

The  actual  feasible  range  of  a  nuclear  weapons  problem 
would  include  these  harder  targets.  Therefore,  weapon  one 
is  a  superior  weapon  to  weapon  two  in  a  realistic  scenario. 

Returning  to  economics,  if  weapon  two  is  really  an 
inferior  good,  the  partial  of  W2  with  TC  (equation  4.15a) 
should  be  negative. 


*2 

bTC 


2#>11C2  *12C1 


2(*22C1  +  *11C2  “  ^12C2C1> 


(4.  15a) 


Setting  C2  equal  to  one  and  substituting  for  the  p 
values,  one  gets  equation  4.15b.  Setting  C2  equal  to  one 
does  not  change  the  problem  in  any  way  since  it  is  only 
changing  the  units  of  measure. 


bW2  -2.8  +  Z.9Cl 

bTC  -7.6C2  +  7-BCj  -  2.8 


(4. 15b) 


Since  weapon  one  is  a  better  weapon  than  weapon  two 
over  most  of  the  range,  the  only  feasible  solutions  will  be 
where  the  price  of  weapon  two  is  less  than  the  price  of 
weapon  one.  Since  the  price  of  weapon  one  is  greater  than 
weapon  two,  and  the  price  of  weapon  two  is  one,  the  numera¬ 
tor  of  equation  4.15b  is  positive.  Now  in  order  for  equa¬ 


tion  4.15b  to  be  negative,  the  denominator  must  be  negative 


It  can  be  easily  seen  that  -for  all  values  o-f  greater  than 
one  (the  price  o-f  C2> ,  the  denominator  is  negative.  Hence, 
the  partial  (equation  4. 15b>  is  negative  and  weapon  two  is 
an  inferior  good. 

How  does  this  apply  to  the  weapon  problem?  It  has 
already  been  shown  that  weapon  one  is  a  better  weapon  com¬ 
pared  to  weapon  two.  Yet,  i-f  weapon  one  costs  more  than 
weapon  two,  not  as  many  can  be  bought  as  weapon  two.  At 
lower  -funding  levels,  the  most  payoff  can  be  accomplished 
using  a  larger  number  of  cheaper  weapons.  As  the  funding 
level  increases,  more  weapons  will  be  exhausting  the  tar¬ 
gets,  and  weapon  one  will  give  a  better  payoff.  In  other 
words,  at  lower  funding  levels  quantity  is  more  important 
than  quality,  and  at  higher  levels  quality  is  more 
important. 

As  for  the  effect  of  reducing  the  cost  of  a  weapon, 
table  4.5  illustrates  this  with  an  example. 


FUNDING:  4000 

Cl* 

5.0 

C  : 

4.0 

BD— POL YNOM I AL 

wj: 

244 

w|: 

695 

VD— POLYNOMIAL 

WJ: 

388 

515 

FUNDING:  4000 

Cl* 

5.0 

C2* 

3.95 

BD— POL YNOM I AL 

wj: 

207 

W?: 

750 

VD— POL YNOM I AL 

W* 

*v 

380 

4 

550 

FUNDING:  4000 

C.: 

4.9 

C-S 

4.0 

BO-POLYNOMIAL 

wj: 

368 

W?: 

549 

VD— POL YNOM I AL 

w}, 

427 

W  : 
*2* 

476 

TABLE  4.5 


RATE  OF  TECHNICAL  SUBSTITUTION  OF  W  FOR  W.  AT  MEAN  (RTS  . ) 

j  »  Ji 

rts21  rts31  rts32 

Polynomial-BD  1.310  2.471  1.887 

INTENSITY  OF  I  DIVIDED  BY  J  (I  .) 

I  I  i  -i1 

21  *31  32 

CES  1.309  2.492  1.904 

TABLE  5. 3 

the  production  process.  Since  there  is  perfect  ease  of 
substitution,  the  inputs  can  be  freely  substituted  for  each 
other.  If  the  inputs  can  be  freely  substituted,  the  tech¬ 
nology  of  the  production  process  will  take  the  better  pro¬ 
ducing  input,  hence,  the  intensities  measure  comparable 
marginal  output  (or  product).  If  the  inputs  were  all  priced 
the  same,  one  would  only  use  the  best  output  produ-  nw 
input.  This  all  hinges  on  the  fact  that  a  mix  of  inputs  is 
not  needed  for  production,  i.e. ,  the  perfect  ease  of  substi¬ 
tution  assumption. 

A  very  interesting  implication  is  that  the  intensity 
parameters  of  the  CES  can  now  be  used  to  make  decisions 
regarding  the  mix  of  weapons  to  buy.  The  critical  points 
are  given  in  equations  5.3a-5.3c.  Note  that  C^  refers  to 
the  price  of  weapon  i,  and  the  a,  p,  and  l-«-p  terms  are  the 
intensity  parameters  of  the  CES. 

P 

- C. 


C„  = 


(5.3a) 


(5.2) 


Y  *  17 ( . 46Wj  +  .35W2  +  . 19W3> - 

The  most  interesting  observation  is  that  the  returns  to 
scale  is  much  less  than  one.  This  problem  is  characterized 
by  noticeably  diminishing  returns  to  scale.  This  agrees 
Mith  earlier  observations  in  the  chapter  that  there  are 
fewer  targets  for  the  total  weapons.  Many  of  the  weapons  in 
this  problem  are  attacking  very  difficult  targets  or  are 
being  doubled  up  on  other  targets.  The  effect  is  that  there 
is  considerably  less  return  at  higher  levels  of  output.  The 
CES  assumes  that  the  returns  to  scale  is  constant  at  0.65 
across  the  range  of  input  combinations.  This  is  not  really 
true,  the  CES  has  in  effect  averaged  the  returns  to  scale 
across  the  range. 

The  distribution  parameters,  or  intensity  parameters, 
have  the  same  interpretation  as  in  the  two  variable  case. 
They  measure  comparative  mean  marginal  products.  The  term 
comparative  is  used  because  they  are  not  the  mean  marginal 
products  of  one  weapon,  but  of  some  unknown  units  of  a 
weapon.  Dividing  the  intensities  by  each  other  is  the  same 
as  computing  the  marginal  products  and  dividing  them  by  each 
other.  The  result  is  the  rate  of  technical  substitution. 
Table  5.3  illustrates  this  property. 

The  fact  that  there  is  perfect  ease  of  substitution  is 
why  the  intensities  exhibit  mean  marginal  values.  In  more 
typical  economic  problems,  the  intensity  of  an  input  is  a 
measure  of  how  intensive  that  input  is  in  the  technology  of 


other  near  the  Man,  and  not  more  than  twenty-five  percent 
different  at  the  extreae  points. 


ECONOMIC  VALUES  AT  THE  MEAN 

MARGINAL  PRODUCTS  OUTPUT  ELASTICITIES 


FUNCTION 

u 

W1 

W2 

W 

*3 

W1 

W2 

W3 

Pol ynomi al -BD 

.72 

.55 

.29 

.21 

.25 

.20 

Pol ynomi al -VD 

.70 

.56 

.33 

.21 

.26 

.23 

CES-BD 

.70 

.56 

.28 

.21 

.25 

.20 

TABLE  5.2 

The  next  two  sections  will  sore  closely  evaluate  the 
econoeic  interpretation  of  these  functions. 

1*  IHE  CES  PRODUCTION  FUNCTION-  Since  the  CES  is  the 
easiest  to  interpret  in  economic  terms,  it  is  evaluated 
first.  The  form  of  the  CES  for  three  variables  is  given  in 
equation  5.1 

— p  —p  —p  —  u/p 

Y  -  ACaWj  +  +  <l-«-P>W3  )  (5.1) 

The  interpretation  of  the  CES  parameters  are  the  same 
as  in  the  two  variable  case.  Using  a  nonlinear  estimation 
technique  to  evaluate  the  parameters  resulted  in  the  P 
parameter  taking  on  the  value,  P  *  -1.05,  which  is  outside 
the  range  for  P.  A  P  value  of  negative  one  implies  perfect 
ease  of  substitution.  Therefore,  P  was  set  to  its  lower 
limit,  negative  one,  and  the  nonlinear  estimation  technique 
was  applied  on  the  remaining  parameters.  The  estimated 
equation  is  given  below  (equation  5.2).  The  function  was 
estimated  using  the  bias  design  points. 


5.3 
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•V 
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is  obvious  sines  thsrs  art  three  variables  instead  of  two. 
But  also,  the  number  of  targets  (not  target  types)  was 
increased  for  the  three  variable  case.  Because  the  increase 
in  weapons  was  proportionately  greater  than  the  increase  in 
targets,  the  result  is  that  there  are  less  targets  per  total 
weapons  involved  than  in  the  two  variable  problem.  The 
effect  in  the  model  is  to  increase  some  right  hand  side 
(RHS)  values.  Therefore,  the  two  problems  are  entirely 
different,  although  both  use  the  same  methodology  to  maxi¬ 
mize  damage  expectancy.  This  poses  no  problems  in  the 
analysis,  since  the  objective  is  not  to  analyze  some  speci¬ 
fic  weapons,  but  rather  to  study  this  methodology  of  apply¬ 
ing  RSM  and  economics,  and  the  information  gained  through 
using  this  methodology  to  determine  cost  effective  weapon 
choices. 

B.  INTERPRETATION  VALUE 

A  simple  comparison  at  the  mean  shows  that  the  func¬ 
tions  give  similar  economic  values  (Table  5.2).  A  more 
convincing  argument  is  to  look  at  the  economic  values  over  a 
wide  range  of  the  inputs.  This  was  also  accomplished  during 
the  research.  The  results  over  the  range  of  input  combina¬ 
tions  for  both  functions  agreed  quite  closely  on  these  basic 
economic  measures.  ’Quite  closely’  is  a  subjective  measure 
of  the  agreement  among  the  functions.  A  more  objective 
judgement  is  that  the  values  are  within  five  percent  of  each 
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Similar  to  thm  two  variablm  analysis,  for  the  thrsm 
variable  case,  there  are  two  important  elements  in  choosing 
a  function  to  analyze  the  surface,  prediction  ability  and 
interpretational  ability.  Zn  terms  of  prediction  value,  the 
polynomial  is  clearly  the  better  fitting  function  (table 
5.1).  Regardless  of  prediction  measure  of  merit,  the  poly¬ 
nomial  provides  a  better  fit  to  the  surface.  The  next  best 
fitting  function  is  the  CES.  Surprisingly,  all  other  func¬ 
tions  had  percent  fits  around  90  percent.  Based  on  these 
results  and  those  of  chapter  IV  the  polynomial  and  the  CES 
are  the  only  two  functions  evaluated  in  this  chapter. 

RANDOM  DESIGN  PERCENT  FIT 

DATA  POINTS  TO  SURFACE 


EUMCUEW 

SDS 

Polynomial  <VD> 

5446 

0.98 

Polynomial  (BD) 

6013 

0.98 

CES  (BD) 

45706 

0.93 

Linear  (BD) 

86875 

0.90 

Linear  (VD) 

87818 

0.90 

MNH  (BD) 

128150 

0.91 

Cobb-Dougl as  (BD) 

228665 

0.88 

Table  5.1 


Before  preceding  into  the  interpretation  of  the  func¬ 
tions,  the  question  arisesj  why  doesn’t  the  other  produc¬ 
tion  functions  fit  as  well  as  in  the  two  variable  case?  The 
primary  reason  is  that  the  LP  problems  are  different.  This 


First,  the 


closely  on  some  major  economic  interpretations, 
production  process  has  high  ease  of  substitution.  Secondly, 
there  are  diminishing  returns  to  scale.  Thirdly,  the  analy¬ 
sis  leads  one  to  buy  all  of  one  type  of  weapon  for  most  cost 
relationships.  Fourthly,  weapon  one  is  a  generally  superior 
weapon  compared  to  weapon  two. 

Of  the  production  functions,  if  prediction  and  inter — 
pretation  are  of  similar  importance,  the  CES  would  seem  to 
be  the  best  to  use.  It  fits  very  well,  98%  to  the  random 
surface,  and  has  very  descriptive  properties. 


The  C-D  does  fit  fairly  well  away  from  the  axis.  The 
elasticities  are  good  measures  of  the  mean  percent  change  in 
output  for  one  percent  changes  in  the  input.  Still,  it 
would  appear  that  one  of  the  preceding  functions  could  give 
a  better  fit  and  provide  as  much  economic  insight  into  the 
problem,  if  not  more. 


6.  IHE  MULTIPLICATIVE  N0N-H0H06ENE0US  PRODUCTION 
FUNCTION.  The  MNH  is  similar  to  the  Cobb-Dougl as,  but  also 
allows  for  variable  elasticities  of  substitution.  The  form 
of  the  MNH  is  given  in  equation  4.24. 


Y 


^(a  +  6LN  (Wj  >  ) 


(4.24) 


The  MNH,  when  estimated  using  the  variance  design,  was 
an  extremely  poor  fit  to  the  surface  of  the  random  data 
points;  on  average  the  error  was  45  percent.  The  bias 
design  was  somewhat  better;  the  fit  was  88  percent,  which 
translates  into  a  12  percent  error  on  average.  The  esti- 
mated  parameters  of  the  MNH  for  the  bias  design  are  given  in 
equation  4.25. 


Y  -  0.001wJ1-83-°-26LN(W2)>WJ-9° 
Since  the  function  is  not  a  very  good  fit 
surface,  it  will  not  be  evaluated  any  further. 


(4.25) 
to  the 


C.  SUMMARY 


In  conclusion,  a  number  of  observations  can  be  made. 
Among  the  functions  that  fit  the  surface  well,  the  poly¬ 
nomial,  the  CES,  the  linear,  and  the  ANH,  all  correlate 
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elasticities  are  constant,  independent  o-f  the  production 
level.  The  output  elasticities  added  together  measure  the 
returns  to  scale  <.87).  The  returns  to  scale  in  the  C-D 
case  is  a  little  less  than  the  CES  measured  (.94).  They  do 
both  agree  that  there  is  diminishing  returns  to  scale. 

Figure  4.6  is  a  three  dimensional  plot  of  equation 
4.23.  The  most  noticeable  inconsistency  of  the  C-D  is  seen 
in  figure  4.6,  that  is,  the  C-D  requires  positive  levels  of 
both  inputs  before  there  is  any  output.  Not  only  does  this 
disagree  with  the  problem  but  also  with  the  implications  of 
the  preceding  production  functions.  The  main  implication  of 
the  other  production  functions  has  been  that  there  is  per¬ 
fect  ease  of  substitution  or  at  least  very  high.  The  C-D 
does  not  agree  with  this.  It  implies  that  when  the  costs  are 
varied,  one  always  buys  a  mix  of  weapons.  This  can  be  seen 
from  figure  4.6)  one  never  buys  zero  of  a  weapon. 


It  is  not  necessary  to  set  up  the  Lagrangian  in  the 
sieple  case  o-f  a  linear  eodel.  Since  the  substitution  is 
infinite,  it  is  eerely  a  natter  of  finding  the  slope  of  the 
isoquants  to  coapute  the  ratio  of  costs  needed  to  buy  a 
Mixture  or  buy  all  of  one  weapon.  If  the  ratio  of  the 
prices  C]/C2  aKCMds  tw^ara  ***  represents  the  mar¬ 
ginal  products),  then  one  buys  all  weapon  two.  If  the  ratio 
of  is  less  than  MP^/MP^  one  buys  all  weapon  one.  If 
the  ratios  of  costs  and  marginal  products  are  equal,  then 
any  combination  of  weapons  at  the  costs  are  appropriate  to 
maximize  damage  expectancy. 

5.  IHE  COBB— DOUGLAS  PRODUCTION  FUNCTION.  The  C-D  also 
has  very  descriptive  properties  and  a  simple  form.  It  has 
been  used  extensively  in  economics.  The  form  of  the  C-D  is 
given  in  equation  4.22. 


Y  ®  AW^W^  (4.22) 

The  C-D  was  estimated  using  the  bias  design.  The  C-D 
is  linear  in  the  logs.  Therefore,  since  the  bias  design  is 
orthogonal  for  functions  of  order  one,  the  parameters  are 
not  correlated.  Equation  4.23  is  the  estimated  function. 

(4.23) 

The  values,  0.33  and  0.54,  are  the  output  elasticities 
for  weapon  one  and  two  respectively.  These  measure  the 


Y  »  3.6wJ,33W2*54 


percent  change  in  output  for  a  one  percent  change  in  the 
input.  In  contrast  to  the  other  functions,  the  output 
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(4.19) 


The  linear  function  was  estimated  using  the  bias  de- 
sign,  which  is  orthogonal  for  models  of  order  one,  and  is 
given  asi 

Y  *  25.2  +  0. 82 +  0.75W2  (4.20) 

The  Pj  and  p2  terms  (0.82  and  0.75)  represent  the 
marginal  products  of  weapons  one  and  two  respectively.  The 
marginal  products  have  a  very  obvious  definition  from  the 
equations  the  increase  in  output  for  a  one  unit  increase  in 
the  input.  Additionally,  the  marginal  products  are  constant. 
The  other  production  functions  have  decreasing  marginal 
products,  a  more  intuitive  result.  If  the  ^  term  had  been 
zero,  then  the  linear  function  would  be  homogeneous  of 
degree  one.  Since  ^  is  equal  to  25.2,  the  model  is  quasi-* 

homogeneous  of  degree  one  for  high  levels  of  the  inputs,  but 
clearly  nonhomogeneous  for  lower  levels  of  the  inputs. 
Additionally,  the  isoquants  for  equation  4.20  are  deary 
linear  and  therefore,  the  substitution  is  infinite.  This 
is  the  same  result  as  in  the  CES  model.  The  polynomial  sup¬ 
ports  high  ease  of  substitution,  but  not  perfect  ease  as  the 
linear  and  CES  are  interpreted. 

The  elasticities  are  easily  computed  by  equation  4.21 
as  the  marginal  product  of  weapon  'i*  times  the  level  of 
weapon  *i*  divided  by  the  output  at  that  level  of  the 
inputs. 

EY,W  *  (4.21) 
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There  is  multi- 


0.085W2LN(W1) 
<4.17) 

0.026W2LN(M1> 
<4.  18) 

4.  JHE  LINEAR  PRODUCTION  FUNCTION.  The  linear  model 
is  probably  the  simplest  o-f  all  the  models  to  understand. 
While  the  coefficients  of  a  linear  model  can  be  obscure  due 
to  differences  in  the  magnitude  of  the  inputs,  this  is  not 
the  case  in  this  study.  A  one  unit  increase  in  either  input 
represent  a  one  weapon  increase.  The  form  of  the  linear  is 
given  below  (equation  4.19). 


that  both  have  a  positive  interaction  term, 
col linearity  in  both  design  estimates. 

BIAS  DESIGN  ESTIMATED 
Y  -  7.8  ♦  0.53WJ  +  1.15W2  +  0.079W1LN(W2) 

VARIANCE  DESIGN  ESTIMATED 
Y  »  27.9  +  0.78WJ  +  0.74W2  +  0. 023W1LN <W2) 


(4.16) 
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Y  -  40  +  41Wl  +  *2W2  +  ^12W1LN(W2>  +  p21W2LN(W1) 

The  ANH  is  basically  a  1 inear  model  with  two  interac¬ 
tion  teres.  This  particular  eodel  was  probably  the  aost 
difficult  to  estimate,  especially  for  aore  than  two  vari¬ 
ables.  The  two  interaction  terms  produced  opposite  signs 
when  estiaated  using  either  the  bias  design  or  the  variance 
design.  In  all  other  cases  of  interactions,  the  effect  has 
been  negative.  The  interpretation  of  this  model  with  a 
nonnegative  interaction  term  is  suspect.  The  results  do  not 
agree  with  intuition  in  certain  places.  The  major  diffi¬ 
culty  is  for  the  case  of  a  large  number  of  weapon  two  and  a 
small  number  of  weapon  one;  the  marginal  product  of  weapon 
one  is  negative.  There  should  be  no  negative  marginal 
products  at  this  level  of  the  inputs,  in  fact,  the  other 
production  functions  estimate  the  marginal  product  of  weapon 
one  at  anywhere  from  0.66  to  0.82.  While  the  ANH  strongly 
agrees  with  the  other  functions  in  the  other  ranges  for 
economic  values,  this  is  one  area  where  it  is  extremely 
suspect. 

Because  of  these  problems  with  the  ANH,  it  will  not  be 
discussed  in  terms  of  economic  meaning.  For  completeness, 
the  actual  estimated  function  is  given  in  equation  4.17  and 
a  graph  of  the  function  in  figure  4.S.  The  variance  design 
estiaated  equation  is  given  in  4.18  for  comparison  purposes. 
Note  that  the  linear  coefficients  are  quite  different,  and 
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does  not  seem  to  agree  Mith  strategic  nuclear  military 
judgement.  The  reason  is  that  the  LP  model  has  as  its 
objective,  damage  expectancy.  It  does  not  model  deterrence. 
Risk  Mould  be  a  main  factor  in  the  results  if  it  did,  and 
Mould  probably  lead  to  conclusions  Mhere  a  combination  of 
Meapons  is  the  best  choice.  What  does  this  say  about  simple 
models  of  this  type  used  in  a  strategic  decision  making 
evironment?  They  give  economic  reasons  for  choosing  one 
Meapon  over  another,  but  a  more  complicated  model  may  be 
needed  to  give  sound  advice  in  force  composition.  This  is 
not  to  say  that  the  model  could  not  be  used  to  examine  force 
increases.  Given  that  the  strategic  force  is  already  a  mix 
of  Meapons,  Mhat  is  the  best  economic  choice  for  increasing 
the  capability  of  the  force?  Here  the  results  imply  that 
this  best  choice  Mill  be  to  buy  all  of  one  Meapon.  The  re¬ 
sults  could  also  be  used  for  examining  the  particular  Meap¬ 
ons  in  one  of  the  legs  of  the  triad.  For  example,  the  B-l 
bomber  and  Stealth  bomber  could  be  analyzed  for  best  choice 
as  an  addition  to  the  present  bomber  leg  of  the  triad. 

3-  IH§  ADDITIVE  NON^HOMOGENEOUS  PRQDyCIION  FUNCTION 
(ANH) .  The  ANH  Mas  originally  developed  for  a  very  specific 
telecommunications  problem  Mhere  there  Mas  some  aprior  in¬ 
formation  that  the  problem  did  not  have  constant  elasti¬ 
cities  of  substitution.  The  returns  to  scale  changed  over 
time  and  the  elasticities  of  output  also  changed.  The  form 
of  the  ANH  is  given  in  equation  4.16. 
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As  one  would  expect,  reducing  the  cost  of  either  weapon 
causes  an  increase  in  the  number  of  that  weapon  purchased. 
Further,  a  reduction  in  the  cost  of  any  weapon  also  causes 
less  of  the  other  weapon  to  be  bought.  A  more  rigorous 


proof  is  to  take  the  partial  of  with  respect  to  either 
or  C2>  The  resulting  equation  is  a  very  difficult  equation 
to  understand  (equation  4.15c). 

<^22C1+*11C2~^12C1C2>  3 

-  tTC<2,..C  12C,)-.2C?^1C,C,1 


2(*22C1  +  *11C2  ^12C2C1> 


(4.15c) 


Since  the  denominator  is  squared  and  hence,  always 
positive,  one  need  only  show  that  the  numerator  is  negative 
to  prove  the  assertion  that  a  reduction  in  price  will  cause 
an  increase  in  purchase.  Unfortunately,  multiplying  the 
numerator  out  yielded  a  very  complicated  expression.  Even 
setting  a  weapon  price  to  one  and  working  in  the  feasible 
range  did  not  simplify  the  expression  enough  to  make  any 
conclusions  about  the  sign. 

One  last  comment  in  regard  to  the  polynomial,  unlike 
the  CES,  there  is  a  range  of  relative  cost  ratios  where  one 
would  buy  a  combination  of  the  two  weapons.  Yet,  as  was 


implied  in  the  plots  of  the  isoquants  (figure  4.3),  this 


range  is  small,  lending  credence  to  the  CES  result  that  it 
is  economically  efficient  to  buy  all  of  one  weapon  to  get 
the  most  destruction  for  a  given  funding  level.  This  result 
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If  the  above  equalities  hold,  any  combination  of  weap¬ 
ons  would  be  the  best  economic  buy.  This  is  true  because 
the  ratio  of  the  prices  for  any  two  inputs  are  equal  to  the 
slope  of  the  isoquants  (the  ratio  of  the  intensities)  for 
those  two  inputs.  The  results  are  similar  to  those  attained 
in  the  two  variable  case.  The  best  economic  decision  is  to 
buy  any  combination  of  weapons,  including  the  choice  to  buy 
all  of  one  weapon.  The  case  where  the  above  equalities  hold 
will  be  rare,  more  likely  will  be  buy  decisions  where  none 


of  the  above  equalities  hold.  These  cases  as  well  as  the 
above  case  are  all  summarized  in  figure  5.1. 
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FIGURE  5.1 


If  the  two  conditions  under  a  particular  buy  decision 
hold,  then  that  decision  is  as  good  as  any  other  possible 
buy  decision.  Other  buy  decisions  could  also  hold  and  they 
would  represent  an  equally  best  economic  choice.  Rarely 
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though,  will  more  that  one  of  the  above  decisions  be  true. 
Most  ratios  of  prices  Mill  not  be  equal  to  the  ratio  o-f  the 
intensities.  In  all  these  cases,  only  one  o-f  the  above 
decisions  Mill  have  both  conditions  holding,  and  that  deci¬ 
sion  Mill  be  the  best  economic  choice. 

The  CES  does  have  some  nice  interpretation  values  in 
this  problem.  The  -fit  is  not  as  good  as  the  polynomial  (93% 
versus  98%),  but  it  alloMs  for  much  easier  interpretation. 
Of  course,  one  must  decide  if  the  CES  percent  fit  of  93% 
captures  the  surface  Mel 1  enough  to  be  able  to  make  deci¬ 
sions  using  it. 

2.  IHE  POLYNOMIAL  PRODUCTION  FUNCTION.  The  polynomial 
function  increases  quite  dramatically  in  terms  Mith  the 
addition  of  the  third  variable.  The  larger  (more  variables) 
a  problem  contains,  the  more  difficult  the  polynomial  is  to 
evaluate.  The  estimated  function,  using  the  variance  de¬ 
sign,  is  given  in  equation  5.6.  The  variance  design  poly¬ 
nomial  Mas  used  instead  of  the  bias  design  because  both  fit 
equally  Mell,  and  the  variance  design  has  no  multi col linear¬ 
ity. 

y  «  -40  +  1.2Wj  +  1.1W2  +  0.9W3 

-4.5(10”4)W1W2  -2.6(10_4)W1W3  -3.  9  ( 10“4) 
-3.8(10_4)W^  -3.2(10"4)W2  -3.3(10”4)W^  (5.6) 

The  coefficients  of  the  linear  terms  are  positive  and 
give  the  marginal  products  of  the  Meapons  at  Iom  input 


levels.  They  also  give  the  maximum  possible  output  -for  a 
weapon.  As  the  level  of  the  inputs  increase,  the  interac¬ 
tion  terms  and  quadratic  terms  begin  to  decrease  the  return 
on  a  weapon.  In  fact  at  a  level  of  weapons  beyond  the 
feasible  range  of  weapons  used  in  this  problem,  the  output 
would  actually  begin  to  decrease. 

The  coefficients  of  the  interaction  and  quadratic  terms 
are  also  slightly  affected  by  the  scale  of  the  problem.  The 
feasible  range  of  weapon  one  is  0-450;  weapon  two  is  0-750; 
weapon  three  is  0-1040.  As  can  be  seen  from  the  ranges,  at 
maximum  level  of  the  inputs,  weapon  three  has  over  twice  as 
many  weapons  as  weapon  one.  Therefore,  even  though  the 
coefficients  of  the  higher  order  terms  involving  weapon 
three  are  smaller  than  the  coefficient  of  weapon  one,  at  the 
higher  input  levels  weapon  three  has  a  larger  negative 
effect  on  its  marginal  output  than  weapon  one.  The  same 
argument  can  be  used  in  comparing  weapon  two  to  weapon  one. 
The  bottom  line  is  that  due  to  the  scale  effects  involved  in 
the  coefficients  of  the  polynomial  function,  it  is  difficult 
to  interpret  the  function  without  a  lot  more  work  and  a  lot 
of  examples. 

In  order  to  evaluate  the  results  of  adding  cost  to  the 
polynomial  function,  a  lagrangian  needs  to  be  set  up.  Due 
to  the  complexity  of  the  problem,  a  more  simplified  approach 
will  be  used  to  study  the  effect  of  costs.  By  holding  one 
of  the  variables  at  some  constant  level,  the  problem  trans- 


fers  into  a  two  variable  problem.  The  results  of  Chapter  IV 
can  then  be  used. 

Since  the  CES  confirmed  that  the  problem  exhibits  per — 
feet  ease  of  substitution,  and  intution  of  a  weapons  problem 
suggests  perfect  ease  of  substitution,  this  will  be  assumed. 
This  assumption  implies  that  the  isoquants  are  flat 
(straight  lines),  meaning  that  in  all  cases,  only  one  weapon 
will  be  bought.  The  actual  estimated  values  for  the  elasti¬ 
cities  of  substitution  for  the  polynomial  were  high;  so  this 
agrees  with  the  assumption. 

Once  perfect  ease  (or  at  least  high  ease)  of  substitu¬ 
tion  is  assumed,  most  buy  decisions  will  buy  only  one 
weapon.  Therefore,  in  analyzing  two  weapons  against  each 
other,  the  other  weapon  can  assumed  to  be  zero. 

The  most  interesting  conclusion  from  the  two  variable 
case  is  that  one  of  the  weapons  was  an  inferior  good.  This 
conclusion  is  tested  in  the  context  of  the  three  weapon 
problem  to  see  if  it  also  holds  here. 

The  most  likely  candidate  for  an  inferior  good  is 
weapon  three.  It  is  definitely  a  poorer  weapon  in  capabil¬ 
ity  than  weapons  one  and  two.  Setting  weapon  two  to  zero 
and  evaluating  the  partial  developed  in  the  previous  chapter 
(equation  4.15a)  gives  equation  5.7.  Note  that  is  the 
price  of  weapon  i,  and  TC  is  the  total  funding. 
bW3  —7 . 6C3  +  2.6C1 

fcTC  2(-3.8c|  -  3.3C *  +  2.60^3) 


(5.7) 


Restricting  <  C^,  keeps  the  problem  within  the 


feasible  range.  Since  weapon  one  is  superior  to  weapon 
three,  none  of  weapon  three  would  be  bought  if  it  cost  more 
than  weapon  one.  In  order  for  equation  5.7  to  be  negative 
and  prove  the  inferior  good  hypothesis,  either  the  numerator 
is  positive  and  the  denominator  is  negative  or  vice  a  versa. 
Changing  the  units  of  cost  such  that  Cj  »  t,  the  following 
relationship  must  hold  for  the  inferior  good  hypothesis: 


**************  ************** 

*  -7.6C-  +  2.6  >0  *  *  -7.6C_  +  2.6  <0  * 

*  3  3  * 

*  AND  *  OR  *  AND  * 

*  -3.8CT  +  2. 6C_  -3.3  <  0  *  *  -3.8C^  +  2.6C  -3.3  >  0  * 

*  *  *  t  *  *  *  *  *  *****  ***  t  t  *  *  r*  ***** 

The  first  set  of  inequalities  holds  for  <  2.6/7. 6  or 
0.34.  The  second  set  of  inequalities  never  holds.  The 
conclusion  is  that  weapon  three  is  an  inferior  good  at 
prices  less  than  34  percent  of  weapon  one.  One  question 
stands  out,  why  does  weapon  three  lose  its  inferior  good 
quality  as  its  price  increases?  It  appears  that  at  levels 
where  the  price  is  higher  than  34  percent  of  weapon  one,  the 
problem  is  infeasible,  that  is,  negative  amounts  of  weapons 
are  being  bought.  To  test  this  possibility  is  beyond  the 
scope  of  this  thesis.  One  interesting  point  to  be  included 
here  is  that  the  CES  would  not  buy  weapon  three  at  prices 


greater  than  40  percent  of  weapon  one  (Intensity  of 
divided  by  the  Intensity  of  W^).  There  may  be  some  correla¬ 
tion  between  these  two  results.  The  difference  may  be  that 


the  polynomial  does  not  have  perfect  ease  of  substitution, 
but  only  high  ease  of  substitution. 

Performing  exactly  the  same  above  analysis  for  weapon 
two  and  weapon  one,  one  finds  that  weapon  two  is  an  inferior 
good  for  prices  of  weapon  two  59  percent  or  less  than  weapon 
one.  For  prices  greater  than  this  ratio,  weapon  two  is  not 
an  inferior  good.  Again,  it  could  be  that  prices  greater 
than  this  ratio  are  infeasible  points  in  the  problem;  this 
was  not  determined.  The  corresponding  ratio  of  prices  break 
point  for  the  CES  was  76  percent  (Intensity  of  W2  divided  by 
the  Intensity  of  W^). 

C.  SUMMARY 

The  similarities  between  the  polynomial  and  the  CES  are 
that  they  both  advocate  buying  a  best  weapon  as  opposed  to  a 
mix,  and  that  they  both  support  the  fact  that  the  weapons 
are  very  interchangeable,  substitutable.  Perhaps  this  is  a 
fault  of  the  model.  The  model  may  be  so  simple  that  it  has 
not  captured  the  essential  aspects  of  a  strategic  weapons 
problem,  in  which  case,  the  analysis  has  revealed  a  poten¬ 
tial  limitation  of  using  this  type  of  model. 

There  are  a  lot  of  similarities  in  the  useful  informa¬ 
tion  gained  from  the  two  and  three  variable  cases.  The  next 
chapter  will  attempt  to  succinctly  summarize  the  value  of 
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This  research  rsprsssnts  a  methodology  for  identifying 
cost  effective  sixes  of  not  only  strategic  forces,  but  of 
potentially  any  type  of  weapon  problem.  In  fact,  the 
results  apply  to  any  problem  using  a  deterministic  eodel  to 
generate  data.  The  methodology  is  simply  the  application  of 
economics  and  Response  Surface  Methodology  <RSM>  to 
deterministic  models. 

In  this  thesis,  the  deterministic  model  reflected  a 
nuclear  exchange  problem.  Using  RSM  to  fit  economic 
production  functions  to  the  surface  of  the  nuclear  exchange 
model  and  Lagrangian  techniques  to  examine  the  effect  of 


cost  on  these  economic  functions,  allowed  the  identification 
of  cost  effective  buys  of  weapons.  In  the  nuclear  exchange 
model  used  in  this  study,  the  cost  effective  choices  were 
easily  determined  by  simple  buy  relationships  that  developed 
out  of  one  of  the  economic  production  functions  (the  CES) . 
These  simple  buy  rules  allow  for  the  determination  of  what 
weapon  mix  to  buy  based  on  only  the  relative  prices  of  the 
weapons  and  the  parameters  of  the  CES  production  function. 


B.  SPECIFIC  BESUWIi 

Fiv*  production  functions,  common  in  economics,  were 
used  in  this  research.  They  are  the  Constant  Elasticity  of 
Substitution  function  (Arrow  and  others.  1961),  the  Cobb- 


Douglas,  a  Multiplicative  Non-Homogeneous  (Vi nod,  1972),  an 
Additive  Nan-Hoeogeneous  (Sudit,  1973) ,  and  a  Linear  -func¬ 
tion.  These  -five  functions  and  a  polynoeial  of  order  two 
were  successfully  fitted  to  the  response  surface  of  a  pro¬ 
duction  process  using  RSM.  The  production  process,  as  men¬ 
tioned  above,  is  an  unclassified  nuclear  exchange  model. 
The  input  variables  represent  a  particular  type  of  warhead; 
the  output  or  production  is  damage  expectancy. 

Of  the  six  functions  fitted  to  the  response  surface, 
three  fit  very  well.  These  functions  are  the  polynomial, 
the  CES,  and  the  Linear.  Interpreting  these  functions  gave 
insight  into  the  production  process  that  was  not  apparent 
at  the  beginning  of  the  research.  Additionally,  the 
coefficients  of  the  economic  production  functions  provide 
information  that  is  obscured  by  the  coefficients  of  the 
polynomial.  The  ease  of  interpretation  of  the  economic 
functions,  such  as  the  CES,  provide  for  a  better  vehicle 
for  analysis  and  decision-making  than  the  polynomial. 

The  polynomial  production  function  provided  a  very  good 
fit  in  both  of  the  cases  studied  in  this  research,  a  two 
variable  case  and  a  three  variable  case.  The  good  fit  of 
the  polynomial  may  be  explained  by  the  fact  that  the  poly¬ 
nomial  is  a  portion  of  the  Taylor  expansion  of  any  function. 
While  the  polynomial  has  good  prediction  value,  the  inter¬ 
pretation  takes  more  work  to  understand  than  the  economic 
functions  used  in  the  study.  The  coefficients  of  the  poly- 


nomial  do  not  have  any  direct  relationship  to  the  coeffi¬ 
cients  of  the  economic  production  functions.  The  best  that 
can  be  said  is  that  the  coefficients  of  the  linear  terms  in 
the  polynomial  represent  a  maximum  marginal  product  of  the 
inputs.  Since  all  the  interaction  terms  and  quadratic  terms 
are  negative,  these  serve  to  reduce  the  marginal  product  as 
the  input  levels  increase.  It  is  very  difficult  to  tell 
exactly  what  degree  of  negative  effect  is  caused  by  the 
terms.  The  best  alternative  is  to  compute  the  first  par- 
tials  to  interpret  the  polynomial.  The  first  partials 
represent  the  marginal  products  and  can  be  easily  computed 
for  the  polynomial. 

The  CEU  production  function  also  fits  very  well  in  the 
context  of  this  problem,  although  not  as  well  as  the  poly¬ 
nomial  in  the  three  variable  case.  The  CES  is  an  easy 
function  for  understanding  the  coeffients.  The  returns  to 
scale  is  evident  in  the  function  and  agrees  with  intuition 
that  the  problem  is  one  of  diminishing  returns,  although  the 
CES  only  gives  the  average  of  the  effect  of  returns  to  scale 
over  the  range.  Also,  the  CES  points  out  clearly  the  high 
ease  of  substitution  in  the  problem  which  is  not  entirely 
evident  from  the  polynomial  function.  This  fact  is  probably 
the  major  cause  why  the  analysis  of  the  results  shows  that 
in  all  cases  the  best  choice,  in  economic  terms,  is  to  buy 
only  one  weapon.  Additionally,  the  coefficients  of  the  CES 
gave  good  comparative  results  of  the  average  marginal  value 


of  a  weapon.  Finally,  these  comparative  values  could  be 
used  to  make  quick  and  simple  cost  decisions  regarding 
weapon  buys. 

The  research  showed  that  economic  production  functions 
can  be  easily  fitted  to  the  response  surface  of  a 
deterministic  model  and  that  the  information  gained  from  the 
fitting  of  these  economic  production  functions  can  provide 
insight  into  the  production  process  that  will  be  useful  to  a 
decision  maker.  The  information  concerning  returns  to 
scale,  ease  of  substitution,  and  the  effect  of  cost  on  the 
problem  is  easily  extracted  from  the  economic  production 
functions.  The  CES  illustrated  these  results  the  best. 

Also,  the  research  was  helpful  because  the  use  of 
different  types  of  production  functions  gave  insight  into 
the  problem  which  may  not  have  been  discovered  had  only  one 
type  of  function  been  used.  For  example,  the  use  of  other 
production  functions  highlighted  passible  model  inconsisten¬ 
cies.  The  result  that  the  best  choice  is  to  buy  only  one 
weapon  is  not  an  obvious  fact  for  a  strategic  weapons  prob¬ 
lem  where  emphasis  has  usually  been  placed  on  a  mix  of 
weapons  (e.g.  the  Triad).  Possibly,  the  model  is  too  simple 
and  needs  to  take  into  account  the  value  of  deterrence  and 
risk,  possibly  through  adding  other  constraints  or  goal 
programming.  Additionally,  the  model  need  not  only  be 
applied  to  problems  of  force  composition,  but  also  where  to 
spend  additional  funds  in  the  pursuit  of  a  more  effective 


-force.  In  other  words,  what  weapon  is  the  best  choice  -for 
increasing  the  e-ff ecti veness  of  the  strategic  force?  The 
results  of  this  research  would  allow  for  the  determination 
of  the  most  cost  effective  strategic  weapon  to  be  added  to 
the  existing  arsenal. 

C.  RECOMMEND AT 1 DNS 

This  work  provides  a  foundation  for  the  use  of  economic 
models  in  the  analysis  of  deterministic  models  of  production 
processes.  The  simple  buy  decision  rules  developed  from  the 
CES  production  function  can  be  a  useful  tool  for  the  deci¬ 
sion  maker  in  cases  where  the  inputs  are  substitutable.  An 
interesting  question  would  be  to  examine  the  results  when 
applied  to  a  more  complex  weapon  problem.  Instead  of 
starting  out  with  only  a  few  variables,  all  the  variables 
could  be  used,  and  an  actual  in  depth  analysis  of  the  re¬ 
sults  using  some  of  these  economic  production  functions 
could  be  accomplished.  The  methodology  as  proposed  in  this 
study  would  choose  an  experimental  design  based  on  the  Box 
and  Bhenken  Designs,  run  the  designs  for  the  deterministic 
model,  fit  the  economic  production  functions  used  in  this 
study,  and  apply  a  Lagrangian  or  nonlinear  optimization 
routine  to  analyze  the  effect  of  a  cost  constraint. 

Further  work  in  this  area  would  include  the  development 


of  bias  minimizing  experimental  designs  which  are  also 
orthogonal.  These  would  be  a  tremendous  help  for  analyzing 
the  results  of  RSM  when  applied  to  deterministic  models. 


This  research  has  provided  a  methodology  -for  the  use  o-f 
economics,  RSM,  and  deterministic  modelling.  The  use  of 
this  methodology  has  great  potential  for  understanding 
problems  more  effectively  and  for  giving  the  decision  maker 
useful  information  regarding  cost  effective  decisions. 


Appendix  A 


This  appendix  lists  the  programs  required  to  generate 
and  run  the  Linear  Programming  model.  A  detailed 
description  o-f  the  model  can  be  -found  in  Graney,  1984, 
appendix  A. 

The  -first  program  is  the  control  cards  required  by  the 
linear  programming  package,  "Multi-Purpose  Optimization 
System",  MPOS,  to  execute  the  problem.  The  second  program 
is  a  -fortran  code  program  to  generate  the  data  matrix  Tor 
the  LP  program. 


tilt  First  Proqraa  tttt 

(THIS  IS  THE  CONTROL  CARD  DECK  FOR  HPOS  AEH 
TITLE  NCN 

ARSENAL  EXHANGE  PROBLEM 
(UNIT 
VARIABLES 
XI  TO  X100 
PACKED 
MAXIMIZE 
CONSTRAINTS  IS 
♦+++♦+♦♦♦+♦♦♦♦♦ 

FORMAT 

(2X, I3,2X, 13, 2X, F10.5) 

REHIND 
READ  HN 
OPTIMIZE 


tilt  Second  Proqran  ttttt 
PROGRAM  CREATE2 

DIMENSION  AIJ (0: 15,0: 100) ,CFI J <0: 15, 0: 100) , CVI J (0: 15, 0: 100) 

C 

C 

C  THIS  PR06RAM  CREATES  THE  DATA  MATRIX  FOR  AN  AGGREGATE  ARSENAL 
C  EXCHANGE  MODEL  TO  BE  SOLVED  BY  ’HPOS’.  THE  MODEL  HAS  TEN 

C  TARGET  CATEGORIES  AW)  FIVE  NEAPON  CLASSES.  THE  MODEL  ALLOCATES 

C  A  SINGLE  TYPE  NEAPON  AGAINST  A  PARTICULAR  TARGET.  THE  CURRENT 


C  OBJECTIVE (RON  0)  HAX1H1ZES  THE  DAHA6E  EXPECTED (DE)  TO  THE 
C  TOTAL  TARGET  BASE.  THE  HODEL  ALLOHS  A  SIN6LE  WEAPON  OR  A 

C  PAIR  OF  SAID  WEAPONS  TO  BE  ALLOCATED  TO  A  S1H6LE  TAROET. 

C 

C  NOTE:  THE  RIGHT-HAND  SIDE  VALUES (RHS)  ARE  IN  COL.  ZERO  AND 

C  THE  OBJECTIVE  COEFF.  ARE  IN  RON  ZERO 

C 

C 

C 

C 

C 

C  INITIALIZE  MATRIX  AIJ  TO  ALL  ZEROS. 

C 

DO  200,  J=0,100 
DO  100,  1=0, IS 
AIJ(I,J)=  0.0 
100  CONTINUE 

200  CONTINUE 
C 
C 

C  BUILDING  TARGET  CONSTRAINT  SET (RONS  1  TO  10  ) 

C 

DO  400,  11=0,80,20 
DO  300,  1=1,10 
AIJII, (211+11)1=  1.0 
AIJII, (211-1+111)=  1.0 
300  CONTINUE 

400  CONTINUE 
C 
C 

C  BUILD  HEAPON  CONSTRAINT  SET(ROHS  11  TO  IS) 

C 

C 

JAM 

JBM9 

C 

DO  800,  1=11,15 
DO  600,  J=JA,JB,2 
AIJ(1,J)=  1.0 
AIJ(I,(J+1)  1=2.0 
600  CONTINUE 
C 

JA=JA+20 

JB=JB+20 

C 

BOO  CONTINUE 
C 
C 

C  SET  MS  (COL.  0  RONS  0  THRU  10)  FOR  OBJECTIVE  AND  NUMBER  OF 
C  TARGETS  IN  EACH  CATEGORY. 


AD-A151  773 


UNCLASSIFIED 


A  METHODOLOGY  FOR  IDENTIFVING  COST  EFFECTIVE  STRATEGIC  2/ 
FORCE  NIXES(U)  AIR  FORCE  INST  OF  TECH  HR I GHT -P AT  TERSON 
AFB  OH  SCHOOL  OF  ENGINEERING  T  U  NANACAPILLI  DEC  84 
AFIT/G0R/0S/84D-8  F/G  15/7  NL 


END 

MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-  1963-A 


0.83678 

AIJ (0,88)= 

0.19496 

AIJ(0,76)= 

0.97336 

AIJ(0,89)= 

0.26052 

AIJ(0,77)= 

0.17125 

AIJ (0,90)* 

0.45317 

AIJ (0,78)* 

0.31317 

AIJ(0,91)= 

AIJ(0,79)= 

AIJ (0,92)= 

0.96008 

AIJ (0,80)= 

0.64099 

AIJ (0,93)= 

AIJ(0,B1)= 

0.72773 

A2J(0,94)= 

0.97750 

AIJ (0,82) = 

0.92587 

AIJ (0,95)* 

AIJ (0,83)= 

AIJ (0,96)= 

0.97730 

AIJ  (0,80* 

0.97750 

AIJ (0,97)* 

AIJ(0,BS>* 

0.48420 

AIJ(0,98)= 

0.84685 

AIJ (0,86)= 

0.73395 

AIJ(0,99I* 

0.78184 

AIJ (0,87)= 

0.10270 

AIJ (0,100)= 

0.95241 

RHS  VALUES  RE8UESTED  FOR  EACH  WEAPON  CATEBORY. 

ANY  COEFFICIENTS  HAY  BE  CHAN6ED  IN  THIS  SECTION.  THE  SESSION 
IS  TERHINATED  WEN  THREE  ZEROS  ARE  ENTERED. 

900  PRINTI, 'INPUT  RON  NUHBER,COLUHN  NUMBER, COEFFICIENT: (I, J.CDEFF)’ 
PRINTS, 'FOR  EACH  NEAPON  RHS  AND  ANY  OTHER  COEFFICIENT  TO  BE’ 
PRINTI, ’CHANGED.  ENTER  ONLY  ONE  SET  AT  A  TINE.  TO  TERMINATE’ 
PRINT!, ’ENTER  THREE  ZEROS  (0  0  0)’ 

READS,  I,J,C 

IF(I . EQ. 0. AND. J.EQ. 0. AND. C.EQ.O)  60  TO  902 

AIJ(I,J)*C 

GOTO  900 

% 

902  PRINT!, ’ENTER  THE  RHS  NEAPON  ALLOCATION  VALUES' 

DO  <6,1=11,13 
READI,C 
AI J  (1 , 0)=C 
<6  CONTINUE 
C 

930  PRINTI, ’THE  FILE  CONTAINING  YOUR  ARSENAL  MODEL  MATRIX' 

PRINTI, ’IS  CALLED-HN  FOR  ALL  FORCES  AND  CV/CF  FOR  OPTIONS.’ 

C 

C 

C 

C  ARSENAL  MATRIX  NRITEN  TO  FILE  NTX.  ARSENAL  MATRIX  ALSO 
C  COPIED  FOR  CV  AND  CF  OPTIONS. 

C 

0PEN(12,FILEa’HN’) 

REMIND  12 
DO  954,  1=0,15 
DO  932,  J=0,100 

IF(AIJ(I,J).NE.O.O)  WRITE ( 12, 951 )  I,J,AIJtI,J) 

931  F0RHATI2X,I3,2X,I3,2X,F10.5> 

CFIJ<I,J)*AIJ(I,J) 

CVIJ(I,J)*AIJ(I,J) 


932  CONTINUE 


934  CONTINUE 
CLOSE (12) 

C 

C  COUNTER-FORCE  MATRIX:  SET  PK*0.0  FOR  CV  TAR6ETS  IN  OBJECTIVE 
C 

QPEN(U,FILES’CF’) 

REMIND  11 
DO  13,  11=0,80,20 
DO  11,  1=1,4 
CFIJ(0,I+II>=  0.0 

11  CONTINUE 

DO  12,  1=17,20 
CFIJ(0,I+U)=  0.0 

12  CONTINUE 

13  CONTINUE 
C 

C 

C 

C  COUNTER-VALUE  MATRIX:  SET  PK=0.0  FOR  CF  TAR6ETS  IN  OBJECTIVE 

C 

0PEN(13,FILE*’CV’) 

REMIND  13 
DO  16,11=0,80,20 
DO  15,1=3,16 
CVIJ(0,I+I1)=  0.0 

15  CONTINUE 

16  CONTINUE 
C 

C 

C 

C 

C 

C 

DO  21,  1=0,13 
DO  20, J=0, 100 

IF (CFI J (I , J) . ME. 0. 0)  NRITE<11,29)  Z,J,CFIJ(I,J) 
IF(CVIJ(I,J).NE.O.O)  WRITE! 13, 29)  1, J.CVIJ (I, J) 

20  CONTINUE 

21  CONTINUE 

29  F0RHAT(2X,I3,2X,13,2X,F10.S) 

C 

CLOSE (11! 

CLOSE (13) 

C 

STOP 

END 


A.S 


deosnglx  § 


This  ssction  contains  ths  experimental  designs  used  in 
the  study.  The  designs  are  explained  in  Chapter  III.  The 
random  data  points  are  also  contained  in  this  appendix.  The 
random  data  points  were  used  to  test  the  -fit  o-f  the  sur-face 
to  a  random  set  of  points  to  the  surface. 


EXPERIMENTAL 


BIAS  DESIGN  (2  VARIABLES) 

W1  W2  RESPONSE 
113  188  255.789550 
113  563  548.390880 
338  188  446.673460 
338  563  718.891660 
66  375  368.951800 
384  375  636.373700 
225  640  686.068030 
225  HO  284.720310 
225  225  381.244030 


BIAS  DESIGN  (3  VARIABLES) 

Ml  W2  M3  RESPONSE 
326  543  753  936.892380 
326  543  287  883.645210 
326  207  753  820.779310 
326  207  287  669.880260 
124  543  753  837.946050 
124  543  287  738.305280 
124  207  753  665.235910 
124  207  287  499.332690 
394  375  520  893.147910 
56  375  520  651.915660 
225  657  520  894.097370 
225  93  520  604.202090 

225  375  911  859.986050 
225  375  129  610.025140 
225  375  520  786.167230 


DESIGN  POINTS 


VARIANCE  DESIGN  (2  VARIABLES) 

Ml  M2  RESPONSE 
000  750  574.755680 
450  750  B24. 277880 
000  OOO  0.000000 
450  OOO  382.358800 
000  375  312.936640 
450  375  688.913240 
225  000  191.243410 
225  750  718.455880 
225  375  503.584490 


VARIANCE  DESIGN  (3  VARIABLES) 

Ml  M2  M3  RESPONSE 
450  750  520  1023.035450 
450  000  520  723.572870 

000  OOO  520  345.113760 

000  750  520  797.826350 

000  375  000  314.400160 

000  375  1040  723.161960 
450  375  1040  981.813950 
450  375  000  695.335010 

225  000  000  191.250000 

225  750  000  787.683080 

225  750  1040  965.302085 
225  000  1040  669.632480 

225  375  520  786.167230 

224  374  519  784.606740 

226  376  521  787.634740 


RANDOM  DATA  POINTS 


TWO  VARIABLES  THREE  VARIABLES 


W1 

W2 

RESPONSE 

W1 

W2 

W3 

RESPONSE 

74 

614 

551.81298 

47 

113 

16 

149.59712 

84 

615 

560.84159 

100 

349 

265 

575.78645 

329 

263 

500.29153 

108 

362 

234 

578.78938 

259 

125 

326.36002 

189 

698 

64 

766. 14644 

137 

630 

615.51921 

169 

299 

851 

782.93449 

75 

263 

284.97425 

351 

52 

114 

439.02531 

435 

152 

498.35396 

448 

547 

587 

978.43160 

175 

481 

542.30104 

433 

689 

56 

917.37715 

142 

573 

579.57019 

403 

107 

661 

795.06902 

355 

146 

425.54636 

384 

276 

450 

821.86015 

18 

250 

225.89286 

130 

521 

917 

855.74153 

335 

250 

494.76924 

286 

307 

501 

785.47078 

41 

7 

40.8 

42 

704 

547 

812.63913 

190 

93 

240.53123 

46 

458 

910 

781.90820 

73 

60 

113.05 

32 

730 

738 

838.28725 

230 

96 

539 

617.68826 

11 

160 

545 

481.68624 

5 

406 

347 

570. 93085 

235 

404 

482 

799.21689 

32 

732 

351 

785. 572485 

219 

684 

893 

934. 53033 

244 

439 

233 

746.9563 

147 

243 

58 

378.72357 

132 

202 

911 

707.60477 

11 

248 

300 

441.56219 

Appendix  C 


mtimtmstmmtmtmumttmttmt 
I  BIAS  DESIGN  ESTIMATION  (  TNG  VARIABLES  )  I 
tttmmmuumttmmmttmmmm 


ASD  COMPUTER  CENTER 
NRI6HT-PATTERS0N  AFB.OHIO 

SPSS--  STATISTICAL  PACKAGE  FOR  THE  SOCIAL  SCIENCES 
VERSION  9.0  (NOS)  -  MARCH  06,  1984 


376S00  CM  MAXIMUM  FIELD  LENGTH  REQUEST 


RUN  NAME 

RESPONSE  SURFACE  FIT 

VARIABLE  LIST 

N1,N2,AF 

INPUT  FORMAT 

FREEF1ELD 

INPUT  MEDIUM 

CARD 

N  OF  CASES 

UNKNONN 

COMPUTE 

N113M1IM1 

COMPUTE 

N22=N2tN2 

COMPUTE 

N12=N1IN2 

COMPUTE 

AFL=LN(AF) 

COMPUTE 

N1L*LN(N1) 

COMPUTE 

N2L=LN(N2) 

COMPUTE 

H1L2L=LN (Ml ) SLN (M2) 

COMPUTE 

N12L*N1ILN(N2) 

COMPUTE 

N21L=N2ILN(N1) 

CPU  TIME  REQUIRED..  .051  SECONDS 


RE6RESSI0N  VARIABLES«AF,N1,M2,M12,N11,N22, 

AFL, NIL, W2L, U1L2L, N12L, W21L/ 
RE6RESSIDN  *  AF  WITH  N1.W2/ 

RE6RESSI0N  ■  AFL  NIDI  H1L.H2L, H1L2L/ 
REGRESSION  »  AFL  KITH  H1L,H2L/ 
REGRESSION  *  AF  WITH  N1,N2,N12L,N21L/ 
STATISTICS  ALL 

READ  INPUT  DATA 


VARIABLE 

NEAN 

STANDARD  DEV 

AF 

480.7893 

173.2979 

HI 

225.2222 

112.4653 

N2 

358.5556 

194.0175 

N12 

80758.4444 

62560.9798 

Nil 

61968.1111 

51623.6080 

H22 

162022.3333 

148832.2058 

H12L  1290.5691 

H21L  1889.6317 

CORRELATION  COEFFICIENTS. 


674.2853 

1059.8697 


A  VALUE  OF  99.00000  IS  PRINTER 
IF  A  COEFFICIENT  CANNOT  BE  CONFUTED. 

HI  .53360 

N2  .84472  .00019 

N12  .94504  .67472  .69864 

Nil  .54738  .98150  .02127  .67695 

H22  .82662  .00021  .98561  .68859  -.00575 

AFL  .99075  .52688  .83340  .90573  .54571  .80579 

NIL  .49373  .97286  -.02241  .64061  .91244  .00955  .47910 

N2L  .83137  .00020  .97564  .68159  .05060  .92651  .83318  -.05680 

N1L2L  .95225  .73308  .66508  .96986  .72288  .65487  .93701  .71418  .65668 

N12L  .68963  .97675  .19929  .81351  .96895  .18928  .67577  .93936  .20424  .85360 

N21L  .93108  .20738  .97464  .83701  .21541  .96743  .90900  .19072  .94403  .B06B9  .39978 

AF  HI  H2  N12  Nil  N22  AFL  NIL  N2L  N1L2L  N12L 


1  RESPONSE  SURFACE  FIT 


84/11/21.  16.44.31.  PA6E  6 


Otltttttttttttttf  MULTIPLE  REGRESSION  ttttttttttttttt 
ODEPENDENT  VARIABLE..  AF 


KEAN  RESPONSE 


480.78927  STD.  DEV. 


173.29790 


VARIABLE(S)  ENTERED  ON  STEP  NUNBER  1..  HI 

N2 


MULTIPLE  R  .99906 
R  S9UARE  .99811 
ADJUSTED  R  SQUARE  .99748 
STD  DEVIATION  8.69151 


ANALYSIS  OF  VARIANCE  DF  SUN  OF  SQUARES 


REGRESSION 

RESIDUAL 

COEFF  OF  VARIABILITY 


2.  239804.02954 

6.  453.25450 

1.8  PCT  F 

1587. 21445 


NEAN  SQUARE 
119902.01477 
75.54242 
SIGNIFICANCE 
.000 


VARIABLES  IN  THE  EQUATION  — 


VARIABLE  B  STD  ERROR  B  F  BETA 

SIGNIFICANCE  ELASTICITY 

HI  .82197447  . 27323221E-01  903.00895  .5334374 

.000  .38505 

N2  .75442105  . 15838338E-01  2268.8635  .8444200 

.000  .54242 

(CONSTANT)  25.140493  8.8599901  8.0443953 

.030 


ALL  VARIABLES  ARE  IN  THE  EQUATION. 

VARIANCE/COVARIANCE  MATRIX 

HI  .00075 

N2  -.00000  .00025 

HI  N2 

FILE  NONANE  (CREATION  DATE  >  84/11/21.) 

Otttttttttttttttt  MULTIPLE  RESRESSION  ttttfttttttttt 
ODEPENDENT  VARIABLE..  AFL 

MEAN  RESPONSE  4.11370  STD.  DEV.  .37944 
VARIABLE(S)  ENTERED  ON  STEP  NUMBER  1..  NIL 

N2L 

N1L2L 


MULTIPLE  R 

.99710 

ANALYSIS  OF  VARIANCE 

DF 

SUM  OF  SQUARES 

MEAN  SQUARE 

R  SQUARE 

.99422 

RESRESSION 

3. 

1.14528 

.38174 

ADJUSTED  R  SQUARE 

.99075 

RESIDUAL 

5. 

.00444 

.00133 

STD  DEVIATION 

.03450 

COEFF  OF  VARIABILITY 

.4  PCT 

F 

SIBNIFICANCE 

284.52212  .000 


-  VARIABLES  IN  THE  EQUATION  - 

VARIABLE  B  STD  ERROR  B  F  BETA 

SIGNIFICANCE  ELASTICITY 

NIL  1.8290258  .34358425  28.338259  2.8993432 

.003  1.57B40 

H2L  1.9025979  .31394508  34.722423  3.0432875 

.002  1.78324 

H1L2L  -.25504794  .58471440E-01  19.029232  -3.1452410 

.007  -1.24048 

1.8454840  13.304943 


(CONSTANT)  -4.7320874 


VARIANCE/COVARIANCE  MATRIX  OF  THE  UNNORHALIZED  RE6RESSI0N  COEFFICIENTS, 


NIL  .11805 
N2L  .10744 
N1L2L  -.02005 

.09857 

-.01832  .00342 

NIL 

N2L  N1L2L 

1 RESPONSE  SURFACE  FIT 

84/11/21.  16.44.31.  PAGE 

FILE  NONANE  (CREATION  DATE  >  84/11/21.) 

Ot  1  1  1  1  1  1  1  1  1  1  1  1  I  1  MULTIPLE 
ODEPENDENT  VARIABLE..  AFL 

REGRESSION  ttttttlttttti 

MEAN  RESPONSE 

6.11370  STD.  DEV. 

.37946 

VARIABLE (S)  ENTERED  ON  STEP  NUMBER  1..  NIL 

N2L 


MULTIPLE  R 

.98601 

R  SQUARE 

.97221 

ADJUSTED  R  SQUARE 

.96294 

STD  DEVIATION 

.07305 

ANALYSIS  OF  VARIANCE 

REERESSION 

RESIDUAL 

COEFF  OF  VARIABILITY 


VARIABLES  IN  THE  EQUATION  - 


DF  SUM  OF  SQUARES 
2.  1.11992 
6.  .03202 
1.2  PCT  F 

1W.  93932 


MEAN  SQUARE 
.55996 
.00534 
SIGNIFICANCE 
.000 


VARIABLE 

B 

STD  ERROR  B 

F 

SIGNIFICANCE 

BETA 

ELASTICITY 

NIL 

.33315945 

.43004505E-01 

60.0172B7 

.000 

.52B122B 

.28754 

H2L 

.53611674 

.42340503E-01 

160.32701 

.000 

.8631775 

.50248 

(CONSTANT) 

1.2837147 

.34234712 

14.060594 

.010 

ALL  VARIABLES  ARE  IN  THE  EQUATION. 

VARIANCE/COVARIANCE  MATRIX  OF  THE  UNNORRALIZED  REGRESSION  COEFFICIENTS. 


1 RESPONSE  SURFACE  FIT 


84/11/21.  16.44.31 


PACE  11 


FILE  NONANE  (CREATION  DATE  -  84/11/21.) 

Otlttttttttttttt  MULTIPLE  REGRESSION  tlttttttltlttlt 
ODEPENOENT  VARIABLE..  AF 

MEAN  RESPONSE  480.78927  STD.  DEV.  173.29790 

VARIABLE (S)  ENTERED  ON  STEP  NUMBER  1..  Ml 

N2 

N21L 

N12L 


MULTIPLE  R 

.99981 

ANALYSIS  OF  VARIANCE 

DF 

SUM  OF  SQUARES 

MEAN  SQUARE 

R  SQUARE 

.99963 

REGRESSION 

4. 

240168.16374 

60042.04093 

ADJUSTED  R  SQUARE 

.99926 

RESIDUAL 

4. 

89.12030 

22.28007 

STD  DEVIATION 

4.72018 

COEFF  OF  VARIABILITY 

1.0  PCT  F 

SIGNIFICANCE 

2694.87615  .000 


—  VARIABLES  IN  THE  EQUATION  - 


VARIABLE 

B 

STD  ERROR  B 

F 

SIGNIFICANCE 

BETA 

ELASTICITY 

HI 

.52838286 

.19150545 

7.6126375 

.051 

.3429051 

.24752 

H2 

1.1537230 

.10367612 

123.83563 

.000 

1.2916626 

.86041 

H21L 

-.85208558E-01 

.21158623E-01 

16.217788 

.016 

-.5211256 

-.33489 

H12L 

. 785474 12E-01 

. 35353191E-01 

4.9363567 

.090 

.3056203 

.21084 

(CONSTANT) 

7.7538571 

11.089156 

.48892087 

.523 

ALL  VARIABLES  ARE  IN  THE  EQUATION. 


VARIANCE/COVARIANCE  MATRIX  OF  THE  UNNORHALIZED  RE5RESSI0N  COEFFICIENTS. 


HI 

.03667 

H2 

-.00182 

.01075 

H12L 

-.00663 

.00102 

.00125 

N21L 

.00120 

-.00214 

-.00035  .00045 

HI  H2  H12L  N21L 


C.5 


STEPNISE  ESTIMATION  OF  POLYNOMIAL 


REGRESSION 

STATISTICS 


METHOD«STEPNISE/VARIABLES*AF,  Hl,N2,N12,Nll,N22/ 
REGRESSION  =  AF  WITH  Nl,N2,N12,Nll,N22/ 

ALL 


I  I  *  I  *  I  I  I  I  I  I  I  I  I  I  *  I  *  »  I  I  I  I  I  »  t  *  I  I  $  I  *  I  I  I  I  I  I  *  I  *  I  I  I  *  I  I  I  t  I  I 
VARIABLE (SI  ENTERED  ON  STEP  NUMBER  4..  N22 


MULTIPLE  R 

.99990 

ANALYSIS  OF  VARIANCE 

OF  SUM  OF  SQUARES 

MEAN  SQUARE 

R  SQUARE 

.99979 

REGRESSION 

4.  240207.40101 

60051.83023 

ADJUSTED  R  SQUARE 

.99938 

RESIDUAL 

4.  49.88303 

12.47076 

STD  DEVIATION 

3.33140 

COEFF  OF  VARIABILITY 

.7  PCT  F 

SIGNIFICANCE 

4B15. 41330  .000 


-  VARIABLES  IN  THE  EQUATION  - 


VARIABLE 

B 

STD  ERROR  B 

F 

SIGNIFICANCE 

BETA 

ELASTICITY 

H12 

-.24194505E-03 

. 8370657 4E-04 

8.3343822 

.045 

-.0873428 

-.04064 

H2 

.99273026 

. 42487624E— 01 

345.93085 

.000 

1.1114215 

.74034 

HI 

.91277073 

.33314990E-01 

730.66132 

.000 

.5923617 

.42758 

H22 

-.  2431 1897E-03 

. 49632284E-04 

23.994361 

.008 

-.2087961 

-.08193 

(CONSTANT) 

-21.806124 

9.6083593 

3.1503936 

.086 

VARIABLES  NOT  IN 

THE  EQUATION  ~ 

_  m  m  t _ 

VARIABLE  PARTIAL  TOLERANCE  F 


SIGNIFICANCE 

HU  .16713  .01118  . 8620804  IE-01 

.788 


C.& 


Otttttttttttttttt  MULTIPLE  RE6RESSI0N  tltltlllttlltlt 
09EPEMDENT  VARIABLE. .  AF 

VARIABLE (31  ENTERED  ON  STEP  NUMBER  5..  Nil 


MULTIPLE  R 

.99990 

ANALYSIS  OF  VARIANCE 

DF  SUM  OF  SQUARES 

MEAN  SQUARE 

R  SQUARE 

.99980 

RE6RESSI0N 

5.  240208.79440 

48041.75888 

ADJUSTED  R  SQUARE 

.99946 

RESIDUAL 

3.  48.48963 

16.16321 

STD  DEVIATION 

4.02035 

COEFF  OF  VARIABILITY 

.8  PCT  F 

S18N1FICANCE 

2972.29069 

.000 

-  VARIABLES  IN  THE  EQUATION  - 

VARIABLE 

B 

STD  ERROR  B 

F 

BETA 

SIGNIFICANCE 

ELASTICITY 

H12 

-.241 9237 IE-03 

. 9529653 IE-04 

6.4446997 

-.0873351 

.085 

-.04064 

N2 

.97352646 

.8134B377E-01 

143.21804 

1.0899217 

.001 

.72602 

HI 

.87831110 

.12334072 

50.708855 

.5699984 

.006 

.41144 

H22 

21829319E-03 

. 10169538E-03 

4.6076334 

-.1874752 

.121 

-.07356 

HU 

.7646B259E-04 

.26043978E-03 

.86208041E-01 

.0227791 

.788 

.00986 

(CONSTANT) 

-15.922070 

22.831374 

.48633435 

ALL  VARIABLES  ARE  IN  THE  EQUATION. 

VARIANCE/COVARIANCE  MATRIX  OF  THE  UNNORMALIZED  RE6RESSIDN  COEFFICIENTS. 


HI 

.01521 

H2 

.00844 

.00662 

H12 

-.00000 

-.00000 

.00000 

Nil 

-.00003 

-.00002 

.00000 

.00000 

H22 

-.00001 

-.00001 

.00000 

.00000 

1PA6E  2  BMP  ESTIMATION  OF  CES  PRODUCTION  FUNCTION 


VARIABLES  TO  BE  USED 

1  HI  2  N2  3  H3  4  AF 

OINPUT  FORMAT  IS 
(F3.0,1X,F3.0,1X,F3.0,1X,F10.6) 

(MAXIMUM  LEN6TH  DATA  RECORD  IS  22  CHARACTERS. 

01  N  P  U  T  VARIABLES . 

VARIABLE  RECORD  COLUMNS  FIELD  TYPE  VARIABLE  RECORD  COLUMNS  FIELD  TYPE 
INDEX  NAME  NO.  BE6IN  END  HIDTH  INDEX  NAME  NO.  BEGIN  END  HIDTH 


1  HI 

2  H2 

1PA6E  3  BHDP 


1  1  3  3  F  3  H3 

1  5  7  3  F  4  AF 

ESTIMATION  OF  CES  PRODUCTION  FUNCTION 


1  9  11  3  F 

1  13  22  10.6  F 


REGRESSION  TITLE 

ESTIMATION  OF  CES  PRODUCTION  FUNCTION 


REGRESSION  NUMBER  .  .  0 

DEPENDENT  VARIABLE .  AF 

HEIGHTIN6  VARIABLE . 

NUMBER  OF  PARAMETERS .  4 

NUMBER  OF  CONSTRAINTS  .  0 

TOLERANCE  FOR  PIVOTING . t.OE-08 

TOLERANCE  FOR  CONVERGENCE  .  1.0E-05 

MAXIMUM  NUMBER  OF  ITERATIONS .  SO 

MAXIMUM  NUMBER  OF  INCREMENT  HALVINGS .  S 


PARAMETERS  TO  BE  ESTIMATED 


1  PI 

2  P2  3  P3 

4  P4 

MINIMUM 

-.212676E+38 

-.212676E+38  -.212676E+38 

-.212676E+38 

MAXIMUM 

. 212G76E+38 

.212676E+38  .212676E+38 

.212676E+38 

INITIAL 

12.750000 

.460000  .350000 

.700000 

USIN6  THE  ABOVE  SPECIFICATIONS  THIS  PROGRAM  COULD  USE  UP  TO  450  CASES. 

OBASED  ON  INPUT  FORMAT  SUPPLIED 

1  RECORDS  READ  PER  CASE. 

NUMBER  OF  CASES  READ . 

VARIABLE 

NO.  NAME 

MEAN 

STANDARD 

DEVIATION 

MINIMUM 

MAXIMUM 

1 

HI 

225.000000 

99.545395 

56.000000 

394.000000 

2 

H2 

375.000000 

165.796777 

93.000000 

657.000000 

3 

H3 

520.000000 

229.918308 

129.000000 

911.000000 

4 

AF 

756.770569 

132.705986 

499.332690 

936.892380 

IPAGE  S  WOP  ESTIMATION  OF  CES(P=-1)  PRODUCTION  FUNCTION 


CASE 

NO. 

RESIDUAL 

OBSERVED 

3  AF 

PREDICTED 

3  AF 

COOK 

DISTANCE 

STD.  DEV. 
PREDICTED 

1  HI 

2  H2 

1 

-2.895281 

255.789550 

258.684831 

.021345 

3.809592 

113.000000 

188.000000 

2 

3.716946 

548.390880 

544.673934 

.081861 

4.744735 

113.000000 

563.000000 

3 

.263276 

446.673460 

446.410184 

.000317 

4.467590 

338.000000 

188.000000 

4 

-6.558510 

718.891660 

725.450170 

.305606 

4.931611 

338.000000 

563.000000 

5 

4.990973 

36B. 951800 

363.960827 

.152102 

4.776147 

66.000000 

375.000000 

6 

12.691570 

636.373700 

623.682130 

.745308 

4.480298 

384.000CM 

375.000000 

7 

-5.816733 

686.068030 

691.8B4763 

.193911 

4.709762 

225.000000 

640.000000 

B 

-7.968606 

284.720310 

292.688916 

.211123 

4.113054 

225.000000 

110.000000 

9 

-.027415 

381.244030 

381.271445 

.000001 

3.318322 

225.000000 

225.000000 

imttmttmmmmtttm 

tttl  THREE  VARIABLE  CASE  tilt 

mtmmmmtummmt 

1PA6E  1 

BffDPAR— DERIVATIVE-FREE  NONLINEAR  REGRESSION 
WOP  STATISTICAL  SOFTNARE,  INC. 

1964  HESTNOOD  BLVD.  SUITE  202 
(2131  475-3700 

PROGRAM  REVISED  APRIL  1982 
MANUAL  REVISED  -  1981 

COPYRIGHT  (0  1982  REGENTS  OF  UNIVERSITY  OF  CALIFORNIA 
PROGRAM  CONTROL  INFORMATION 

/PROBLEM  TITLE  IS  ’ESTIMATION  OF  CES  PRODUCTION  FUNCTION’. 
/INPUT  VARIABLES  ARE  4. 

FORMAT  IS  ’ (F3.0,1X,F3.0,1X,F3.0,1I,F10.6)’. 

UNIT  IS  9. 

/VARIABLE  NAMES  ARE  H1,N2,H3,AF. 

/REGRESS  DEPENDENT  IS  AF. 

PARAMETERS  ARE  4. 

/PARAMETER  INITIAL  ME  12.75, .46, .35, .70. 

/END 


BMDP  UNIT  NO.  9  SPECIFIED  IN  THE  INPUT  PARAGRAPH 
HILL  REFER  TO  LOCAL  FILE  NAME  BD5K3C  FOR  THIS  PROBLEM. 


NUMBER  OF  VARIABLES  TO  READ  IN .  4 

TOTAL  NUMBER  OF  VMIABLES .  4 

NUMBER  OF  CASES  TO  READ  IN . TO  END 

INPUT  UNIT  NUMBER  .  9 


REMIND  INPUT  UNIT  PRIM  TO  READING.  .  DATA.  .  .  YES 


VARIABLE  STANDARD 


NO. 

NAME 

MEAN 

DEVIATION 

MINIMUM 

HAHNUR 

1 

HI 

225.222222 

112.465303 

66.000000 

384.000000 

2 

H2 

358.555556 

194.017468 

110.000000 

640.000000 

3 

AF 

480.7B9269 

173.297895 

255.789550 

718.891660 

ITER. 

INCR. 

RESIDUAL  SUM 

PARAMETERS 

NO. 

HALV. 

OF  SQUARES 

1  PI 

2  P2 

3  P3 

USER  ROUTINE  FUN  LOADED  FROM  FILE  L60  AT  ADDRESS  077656B. 

0 

0 

3B0069. 453937 

1.000000 

.550000 

1.000000 

0 

0 

351591.813292 

1.000000 

.500000 

1.000000 

0 

0 

250095.842247 

1.100000 

.500000 

1.000000 

0 

0 

30051.998468 

1.000000 

.500000 

1.100000 

1 

0 

19915.686880 

2.104949 

.532411 

.943083 

2 

0 

0192.384342 

2.247996 

.521871 

.937137 

3 

0 

571.951114 

2.383673 

.521616 

.935671 

4 

0 

458.046663 

2.370562 

.520632 

.939374 

5 

0 

348.665033 

2.337365 

.520569 

.940588 

6 

2 

348.650220 

2.336918 

.520601 

.940628 

7 

0 

348.601419 

2.338077 

.520531 

.940561 

B 

0 

348.601283 

2.338310 

.520535 

.940544 

9 

0 

348.601278 

2.338360 

.520536 

.940540 

10 

2 

348.601278 

2.338361 

.520536 

.940540 

11 

0 

348.601278 

2.338366 

.520536 

.940540 

12 

0 

348.601278 

2.338364 

.520536 

.940540 

1PA6E  4  BNDP  ESTIMATION  OF  CES(P=-1)  PRODUCTION  FUNCTION 


THE  RESIDUAL  SUN  OF  SQUARES  (  *  348.601  )  HAS  SMALLEST  HITH  THE  F0LL0HIN6  PARAMETER  VALUES 


1  PI  2  P2  3  P3 

2.33836  .520536  .940540 


ESTIMATE  OF  ASYMPTOTIC  CORRELATION 

HATRU 

PI 

P2 

P3 

1 

2 

PI 

1 

1.0000 

P2 

2 

-.2113 

1.0000 

P3 

3 

-.9978 

.2518 

1. 

THE  ESTIMATED  MEAN  SQUARE  ERROR  IS  58.10 

ESTIMATES  OF  ASYMPTOTIC  STANDARD  DEVIATIONS  OF  PARAMETER  ESTIMATES  HITH  6  DEBREES  OF  FREEDOM  ARE 

2  P2  3  P3 

8.64469BE-03  1.683183E-02 


1  PI 
227410 


BMP  UNIT  NO.  9  SPECIFIED  IN  THE  INPUT  PARAGRAPH 
MIU  REFER  TO  LOCAL  FILE  NAME  BD4K2C  FOR  THIS  PROBLEM, 


PROBLEM  TITLE  IS 

ESTIMATION  OF  CES(P'-l)  PRODUCTION  FUNCTION 


NUMBER  OF  VARIABLES  TO  READ  IN .  3 

TOTAL  NUMBER  OF  VARIABLES .  3 

NUMBER  OF  CASES  TO  READ  IN . TO  END 

INPUT  UNIT  NUMBER .  9 


REMIND  INPUT  UNIT  PRIOR  TO  READING.  .  DATA.  .  .  YES 
1PA6E  2  BMP  ESTIMATION  OF  CESIP*-1>  PRODUCTION  FUNCTION 


VARIABLES  TO  BE  USED 

1  Ml  2  M2  3  AF 

OINPUT  FORMAT  IS 
(F3.0,1X,F3.0,1X,F10.6) 

OHAXIMUM  LENGTH  DATA  RECORD  IS  IS  CHARACTERS. 
OINPUT  VARIABLES . 


VARIABLE 

RECORD 

COLUMNS 

FIELD 

TYPE 

VARIABLE 

RECORD 

COLUMNS 

FIELD 

INDEX 

NAME 

NO. 

BEGIN  END 

MIDTH 

INDEX  NAME 

NO. 

BEGIN  END 

MIDTH 

1  Ml 

1 

1  3 

3 

F 

3  AF 

1 

9  18 

10.6 

2  M2 

1 

5  7 

3 

F 

1PA6E 

3  BHDP 

ESTIMATION  OF  CES(P=-1>  PRODUCTION  FUNCTION 

REGRESSION  TITLE 

ESTIMATION  OF  CESIP’-l)  PRODUCTION  FUNCTION 


REGRESSION  NUMBER .  0 

DEPENDENT  VARIABLE .  AF 

WEIGHTING  VARIABLE . 

NUMBER  OF  PARAMETERS .  3 

NUMBER  OF  CONSTRAINTS . 0 

TOLERANCE  FOR  PIV0TIN6 . l.OE-OB 

TOLERANCE  FOR  CONVERGENCE . 1.0E-05 

MAXIMUM  NUMER  OF  ITERATIONS .  SO 

MAXIMUM  NUMBER  OF  INCREMENT  HALVINGS .  S 

PARAMETERS  TO  BE  ESTIMATED 

1  PI  2  P2  3  P3 

MINIMUM  -.212676E+38  -.212676E+38  -.212676E+3B 

MAXIMUM  .212676E+38  .212676E+38  .212676E+38 

INITIAL  1.000000  .500000  l. 000000 

USING  THE  ABOVE  SPECIFICATIONS  THIS  PR06RAH  COULD  USE  UP  TO  360  CASES. 

OBASED  ON  INPUT  FORMAT  SUPPLIED  1  RECORDS  READ  PER  CASE. 

NUMBER  OF  CASES  READ .  9 


D.2 


ABBSQdlx  6 


This  section  contains  the  results  o-f  the  CES  production 
function  when  estimated  using  a  nonlinear  estimation 
technique.  The  -First  portion  are  the  results  for  the  two 
variable  case;  the  second  section  are  the  results  of  the 
three  variable  case. 


1PA6E  1 

BHDPAR— DERIVATIVE-FREE  NONLINEAR  REGRESSION 
SHOP  STATISTICAL  SOFTNARE,  INC. 

1964  WESTWOOD  BLVD.  SUITE  202 
1213)  475-5700 

PROGRAH  REVISED  APRIL  1982 
NANUAL  REVISED  --  1981 

COPYRIGHT  (C)  1982  REGENTS  OF  UNIVERSITY  OF  CALIFORNIA 

TO  SEE  REHARKS  AND  A  SUHHARY  OF  NEW  FEATURES  FOR 
THIS  PROGRAH,  STATE  HENS.  IN  THE  PRINT  PARAGRAPH. 


THIS  VERSION  OF  SHOP  HAS  SEEN  CONVERTED  FOR  USE  ON 

CDC  6000  AND  CYBER  SERIES  CONPUTERS  BY 
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PROGRAH  CONTROL  1NF0RHAT1QN 

/PROBLEM  TITLE  IS  ’ESTIMATION  OF  CES(P*-1)  PRODUCTION  FUNCTION’. 
/INPUT  VARIABLES  ARE  3. 

F0RHAT  IS  MF3.0,1I,F3.0,1X,F10.6»’. 

UNIT  IS  9. 

/VARIABLE  NAHES  ARE  W1,W2,AF. 

/REGRESS  DEPENDENT  IS  AF. 

PARAMETERS  ARE  3. 

/PARAHETER  INITIAL  ARE  1,. 5,1. 

/END 


ttmtmmmmmmmmttttmmmmtttt 

t  VARIANCE  DESIGN  ESTIMATION  (  TWEE  VARIABLES  )  t 

ttimtttssittmttstttmtfttmtmtstmmtstt 

Otttttttttttttttt  MULTIPLE  REGRESSION  tttttttttttlttt 
OOEPENDENT  VARIABLE..  A F 

MEAN  RESPONSE  70S. 10239  STD.  DEV.  243.30773 


VARIABLE(S)  ENTERED  ON  STEP  NUMBER  1..  Hi 


MULTIPLE  R 

.99831 

M33 

M2 

N23 

M13 

N12 

N22 

MU 

N3 

ANALYSIS  OF  VARIANCE 

OF 

SUM  OF  SQUARES 

MEAN  SQUARE 

R  SQUARE 

.99662 

REGRESSION 

9. 

825979.88581 

91775.54287 

ADJUSTED  R  SQUARE 

.99054 

RESIDUAL 

S. 

2801.20330 

560.24066 

STD  DEVIATION 

23.66940 

COEFF  OF  VARIABILITY 

3.4  PCT  F 

SIGNIFICANCE 

163.81450  .000 


-  VARIABLES  IN  THE  EQUATION  - 


VARIABLE 

B 

STD  ERROR  B 

F 

SIGNIFICANCE 

BETA 

ELASTICITY 

HI 

1.1681321 

.13749718 

72.176595 

.000 

.8165836 

.37275 

N33 

-.32614660E-03 

.45554629E-04 

51.257883 

.001 

-.5790773 

-.19178 

N2 

1.0879986 

.824980841-01 

173.92778 

.000 

1.2676083 

.57864 

M23 

-.38559420E-03 

.60690778E-04 

40.365987 

.001 

-.3693692 

-.10664 

N13 

-.  26128816E-03 

■  10U5130E-03 

6.6726226 

.049 

-.1501764 

-.04336 

M12 

-.45407410E-03 

. 14026313E-03 

10.480124 

.023 

-.1882074 

-.05434 

M22 

-.31630721E-03 

.87594558E-04 

13.039590 

.015 

-.2920712 

-.09673 

Nil 

-.380638B8E-03 

.24331916E-03 

2.4472237 

.179 

-.1265306 

-.04190 

N3 

.86739792 

.59493754E-01 

212.56610 

.000 

1.4013512 

.63969 

(CONSTANT) 

-39.723798 

27.964356 

2.0178639 

.215 

VARIABLES  IN  THE  EQUATION 


VARIABLE 

B 

STD  ERROR  B 

F 

BETA 

SISNIFICANCE 

ELASTICITY 

H12 

-.60263960E-03 

.32447497E-03 

3.4494744 

-.2466326 

.122 

-.06719 

N3 

.86369212 

.14355706 

36.196280 

1.4963803 

.002 

.59347 

.  ■ 

N2 

1.2460058 

.19905066 

39.184370 

1.5567025 

.002 

.61743 

HI 

1.2405823 

.33206385 

13.957510 

.9305854 

-  •  •  - 

.013 

.368B4 

H23 

-.52343376E-03 

. 14065224E-03 

13.849374 

-.4948856 

.014 

-.13488 

H33 

-.30201 080E-03 

. 11669216E-03 

6.6982537 

-.5555488 

.049 

-.12760 

N22 

— . 40276597E-03 

.22434771E-03 

3.2230146 

-.3852893 

A-  -  1 

.133 

-.08850 

r  .  •« 

N13 

-. 32612592E— 03 

. 23395620E-03 

1.9431284 

-.1850829 

.222 

-.05042 

Nil 

-.30196BB2E-03 

. 62425204E-03 

.23399376 

-.1040510 

.649 

-.02389 

'  •***  > 

(CONSTANT) 

-66.041758 

83.688299 

.62274143 

— 

.466 

»  ... 

ALL  VARIABLES  ARE  IN  THE  EQUATION. 

VARIANCE/COVARIANCE  MATRIX  OF  THE  UNNORMAL I ZED  RESRESSION  COEFFICIENTS. 


HI 

.11027 

H2 

.03813 

.03962 

H3 

.02750 

.01650 

.02061 

N12 

-.00004 

-.00002 

-.00000 

.00000 

. .*•  .* 

N13 

-.00003 

-.00000 

-.00001 

-.00000 

.00000 

__  ^ 

N23 

-.00000 

-.00001 

-.00001 

-.00000 

-.00000 

.00000 

- 

Nil 

-.00018 

-.00006 

-.00005 

.00000 

.00000 

.00000 

.00000 

•  .  \n‘ 

H22 

-.00004 

-.00004 

-.00002 

.00000 

.00000 

.00000 

.00000 

.00000 

H33 

-.00002 

-.00001 

-.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

Ml 

M2 

M3 

M12 

M13 

M23 

Nil 

M22 

M33 

n  -  .-i 

—  VARIABLES  IN  THE  EQUATION  - 


VARIABLE 

B 

STD  ERROR  8 

F 

SISNIFICANCE 

BETA 

ELASTICITY 

HI 

.78473263 

.13327292 

33.632838 

.004 

.5688254 

.37867 

H2 

.74138196 

.B12163B9E-01 

83.374062 

.001 

.8959116 

.39641 

H12L 

.22738666E-01 

. 33438221 E-01 

.46242749 

.534 

.1036654 

.04589 

H21L 

-.261777B1E-01 

.21905433E-01 

1.4281111 

.298 

-.1827965 

-.08088 

(CONSTANT) 

27.938445 

36.914071 

.57282319 

.491 

ALL  VARIABLES  ARE  IN  THE  EQUATION. 

VARIANCE/COVARIANCE  MATRIX  OF  THE  UNNORHALIZED  RE6RESSI0N  COEFFICIENTS. 


HI 

.01830 

H2 

.00597 

.00660 

H12L 

-.00221 

-.00036 

.00112 

H21L 

-.00041 

-.00088 

-.00049  .00048 

HI  N2  N12L  H21L 


mtmmmmtmmmmmmtmmmi 

t  BIAS  DESI6N  ESTIRATION  (  THREE  VARIABLES  )  t 

ttmttmmmmtttumtmitttttstmtm 


Ottlttttttttttftt  MULTIPLE  RE6RESSI0N  ftttttttttttttt 
ODEPENDENT  VARIABLE..  AF 

VARIABLE  (S)  ENTERED  ON  STEP  NUMBER  V..  HU 


MULTIPLE  R  .99734 
R  SQUARE  .99509 
ADJUSTED  R  SQUARE  .98623 
STD  DEVIATION  13.37243 


ANALYSIS  OF  VARIANCE 

DF 

SUM  OF  SQUARES 

MEAN  SQUARE 

RE5RESSI0N 

9. 

245339.80380 

27259.97820 

RESIDUAL 

5. 

1212.50548 

242.50110 

COEFF  OF  VARIABILITY 

2.1  PCT  F 

SISNIFICANCE 

112.41177  .000 


VARIABLES  IN  THE  EQUATION 


VARIABLE 

B 

STD  ERROR  B 

F 

SIGNIFICANCE 

BETA 

ELASTICITY 

NIL 

.37923771 

.17176596 

4.8747104 

.069 

.5213492 

.26571 

N2L 

.41957953 

.15791222 

7.059B769 

.038 

.6274121 

.32004 

(CONSTANT) 

2.2713777 

1.044BS54 

4.7257021 

.073 

ALL  VARIABLES  ARE  IN  THE  EQUATION. 

VARIANCE/COVARIANCE  MATRIX  OF  THE  UNNORMALIZED  REGRESSION  COEFFICIENTS. 


NIL  .02930 

H2L  -.00000  .02494 

NIL  M2L 


1 RESPONSE  SURFACE  FIT 


B4/U/21.  16.45.12. 


FILE  NONANE  (CREATION  DATE  *  84/11/21.) 

Oltttttlltttttttt  MULTIPLE  REGRESSION  ttttltt 
ODEPENDENT  VARIABLE..  AF 

KEAN  RESPONSE  466.28067  STD.  DEV.  268.81649 

VARIABLE (S)  ENTERED  ON  STEP  NUMBER  1..  HI 

N2 

N12L 

N21L 


MULTIPLE  R  .99265 
R  SQUARE  .98535 
ADJUSTED  R  SQUARE  .97071 
STD  DEVIATION  46.00938 


ANALYSIS  OF  VARIANCE 

REGRESSION 

RESIDUAL 

COEFF  OF  VARIABILITY 


DF  SUM  OF  SQUARES 
4.  569630.97839 
4.  B467. 45393 
9.9  PCT  F 


67 


.272 


PAGE  13 


ttttttt 


NEAN  SQUARE 
142407.74460 
2116.86348 
SIGNIFICANCE 
!  .001 


MULTIPLE  R 

.99915 

ANALYSIS  OF  VARIANCE 

DF  SUN  OF  SQUARES 

HEAN  SQUARE 

R  SQUARE 

.99829 

RE6RESS10N 

3. 

35.40401 

11.82134 

ADJUSTED  R  SQUARE 

.99727 

RESIDUAL 

5. 

.00072 

.01214 

STD  DEVIATION 

.11020 

COEFF  OF  VARIABILITY 

2.0  PCT 

F  ! 

973.44024 

SI6N1FICANCE 

.000 

-  VARIABLES  IN  THE  EOUATION  - 


VARIABLE 

B 

STD  ERROR  8 

F 

BETA 

SISNIF1CANCE 

ELASTICITY 

NIL 

.97003215 

•23224050E-01 

1744.0019 

1.3335310 

.000 

.07900 

N2L 

.90220797 

.21339018E-01 

2033.4955 

1.4389132 

.000 

.73398 

N1L2L 

-.14125954 

. 452701 58E-02 

973.00943 

-1.2849180 

.000 

-.41394 

(CONSTANT) 

.10738532E-02 

.10947133 

.23379434E-03 

.988 

ALL  VARIABLES  ARE  IN  THE  EOUATION. 

VARIANCE/COVARIANCE  MATRIX  OF  THE  UM0RMLIZE9  REfiRESSION  COEFFICIENTS. 


NIL 

.00054 

B2L 

.00033 

M64A 

■ NNN IV 

N1L2L 

-.00009 

•  MOM  MM2 

• VN  ■ ***** 

NIL  N2L  N1L2L 


FILE  NONANE  (CREATION  BATE  >  84/11/21.1 

04  4  I  I  I  t  I  I  t  I  I  t  I  t  t  t  t  I  t  t  t  t  1  MULTIPLE  RE6RE88I0N  44444444444444444 
OOCPEHRET  VARIABLE..  AFL 

hear  resume  s.ovsis  std.  kv.  2.10727 

VMIABLE(S)  ENTEREB  ON  STEP  NUHBER  1..  NIL 

N2L 


NULT1FLE  R 

.81575 

ANALYSIS  OF  VARIANT 

jF 

SUH  OF  SQUARES 

NEAR  SQUARE 

r  mm 

.00345 

REfiRESSION 

2. 

23.03990 

11.81998 

ANNBTD  R  MNMN 

.35393 

RESIDUAL 

0. 

11.88477 

1.98079 

STB  BEV1ATI8N 

1.40741 

COEFF  OF  VARIABILITY 

25.7  PCT 

F 

SISNIFICANCE 

3.90729  .037 


1RESP0NSE  SURFACE  FIT  84/11/21.  16.45. 12.  PA6E  7 


FILE  NONANE  (CREATION  DATE  « 
at  ttttttttttttttt 
ODEPENDENT  VARIABLE. .  AF 

84/11/21.) 

MULTIPLE  REGRESSION  ttttttttttttttt 

MEAN  RESPONSE  466.28067 

STD.  DEV.  268.81649 

VARIABLE IS)  ENTERED  ON  STEP  NUMBER  1..  HI 

H2 

MULTIPLE  R  .98997 

R  SQUARE  .98004 

ADJUSTED  R  SQUARE  .97339 

STD  DEVIATION  43.85135 

ANALYSIS  OF  VARIANCE 
REGRESSION 

RESIDUAL 

COEFF  OF  VARIABILITY 

DF  SUN  OF  SQUARES  MEAN  SQUARE 

2.  566560.78681  283280.39340 

6.  11537.64551  1922.94092 

9.4  PCT  F  SIGNIFICANCE 

147.31622  .000 

-  VARIABLES  IN  THE  EQUATION 


VARIABLE 

B 

STD  ERROR  B 

F 

SIGNIFICANCE 

BETA 

ELASTICITY 

HI 

.74656119 

. 79S6SS06E-01 

88.040228 

.000 

.5411562 

.36025 

H2 

.68617210 

.47739303E-01 

206.59222 

.000 

.8289705 

.55184 

(CONSTANT) 

40.989864 

29.234233 

1.9659367 

.210 

ALL  VARIABLES  ARE  IN  THE  EQUATION. 

VARIANCE/COVARIANCE  MATRIX  OF  THE  UNNORMALIZED  RE6RE55I0N  COEFFICIENTS. 


N1  .00633 

N2  0  .00220 

HI  N2 


FILE  NONANE  (CREATION  DATE  *  84/11/21.) 

Otttttttttttlttll  MULTIPLE  REGRESSION  ttttttttttttttt 
ODEPENDENT  VARIABLE..  Aft 

MEAN  RESPONSE  5.48315  STD.  DEV.  2.10727 


VARIABLE (S)  ENTERED  ON  STEP  NUMBER  1..  NIL 

N2L 


1 RESPONSE  SURFACE  FIT 


84/ 11/21.  16.45.12.  PAGE  S 


FILE  NONANE  (CREATION  DATE  *  84/11/21.) 

Otttttttttttttttt  MULTIPLE  REGRESSION  tttttttttlttttt 
ODEPENDENT  VARIABLE..  AF 

MEAN  RESPONSE  466.28067  STD.  DEV.  268.81649 

VARIABLE (S)  ENTERED  ON  STEP  NUMBER  1..  HI 

N2 

H12 

Nil 

N22 


MULTIPLE  R  .99884 
R  SQUARE  .99768 
ADJUSTED  R  SQUARE  .99381 
STD  DEVIATION  21.14243 


ANALYSIS  OF  VARIANCE 

REGRESSION 

RESIDUAL 

COEFF  OF  VARIABILITY 


DF  SUM  OF  SQUARE5 
5.  576757.42575 
3.  1341.00657 
4.5  PCT  F 

258.05575 


MEAN  SQUARE 
115351.48515 
447.00219 
SIGNIFICANCE 
.000 


-  VARIABLES  IN  THE  EQUATION 


VARIABLE 

B 

STD  ERROR  8 

F 

SIGNIFICANCE 

BETA 

ELASTICITY 

HI 

.95834307 

.14607661 

43.040831 

.007 

.6946695 

.46244 

H2 

1.0589761 

.87645966E-01 

145.98508 

.001 

1.2793582 

.85167 

N12 

-.3935B993E-03 

. 12S2B845E-03 

9.8688344 

.052 

-.1747096 

-.07122 

HU 

-. 14263480E-03 

. 2953077 IE-03 

.23329284 

.662 

-.0484257 
-. 02581 

H22 

-. 37B99507E-03 

. 1063107BE-03 

12.709034 

.038 

-.3574225 

-.19050 

(CONSTANT) 

-12.391642 

18.975908 

.42643475 

.560 

VARIANCE/COVARIANCE  MATRIX  OF  THE  UNNORHALIZED  REGRESSION  COEFFICIENTS. 


Ml 

.02134 

M2 

.00132 

.00768 

N12 

-.00001 

-.00000 

Nil 

-.00004 

0 

N22 

0 

-.00001 

Ml  N2 


.00000 
0  .00000 

0  0  .00000 

N12  Nil  N22 


C.9 


mittmtmtmmmmmmtmmmmm 

*  VARIANCE  DES1BN  ESTIMATION  (  TWO  VARIABLES  )  t 

umtmtmmtttttmtmmummmsttm 


1 RESPONSE  SURFACE  FIT 


84/11/21.  16.4S.12.  PA6E  4 


FILE  NONANE  (CREATION  DATE  *  84/11/21.) 

Ottttt*t)tt)ttltt  MULTIPLE  RE6RESS10N  ttltttttttttttt 


VARIABLE 

MEAN 

STANDARD  DEV 

CASES 

AF 

466.2807 

268.8165 

9 

HI 

225.0000 

194.8557 

9 

N2 

375.0000 

324.7595 

9 

N12 

84375.0000 

119324.2693 

9 

Nil 

84375.0000 

91265.5167 

9 

N22 

234375.0000 

253515.3241 

9 

AFL 

5.4831 

2.1073 

9 

NIL 

3.8418 

2.8969 

9 

H2L 

4.1823 

3.1511 

9 

N1L2L 

16.0676 

19.1681 

9 

N12L 

941.0249 

1225.5323 

9 

W21L 

1440.6685 

1877.1157 

9 

CORRELATION  COEFFICIENTS. 

A  VALUE  OF  99.00000  IS  PRINTED 
IF  A  COEFFICIENT  CANNOT  BE  COMPUTED. 


HI 

.54116 

N2 

.82897 

0 

M12 

.79535 

.61237 

.61237 

HU 

.51620 

.96077 

0  .58835 

H22 

.76896 

0 

.96077  .58835  0 

AFL 

.79739 

.48604 

.57848  .41759  .41423  .48822 

NIL 

.49964 

.91317 

.00000  .55920  .76430  .00000  .52135 

H2L 

.79529 

.00000 

.90971  .55708  0  .75886  .62741  .00000 

H1L2L 

.79520 

.57720 

.57454  .89179  .48311  .47926  .46675  .63209 

.63156 

N12L 

.79135 

.66498 

.52629  .94435  .63889  .43901  .43990  .60724 

.57832 

.94287 

H21L 

.80521 

.52848 

.66467  .94640  .44233  .63859  .44228  .57873 

.60466 

.94070 

.88665 

AF 

HI 

H2  H12  HU  H22  AFL  NIL 

M2L 

N1L2L 

N12L 

C.  8 


* 

i 

ITER.  INCR. 

RESIDUAL  SUM 

PARAMETERS 

NO.  HALV. 

OF  SQUARES 

1  PI 

2  P2 

3  P3 

4  P4 

2 

USER  ROUTINE  FUN  LOADED  FROM  FILE  L&O  AT  ADDRESS  Q77656B. 

0  0 

2012146.071836 

12.750000 

.460000 

.350000 

.770000 

0  0 

58356.798237 

14.025000 

.460000 

.350000 

.700000 

0  0 

3242S.418B43 

12.750000 

.506000 

.350000 

.700000 

1 

0  0 

19025.495991 

12.750000 

.460000 

.385000 

.700000 

0  0 

11010.666724 

12.750000 

.460000 

.350000 

.700000 

1  0 

8479.300134 

16.003967 

.462151 

.352606 

.668145 

2  0 

6112.735411 

16.896061 

.462232 

.352563 

.655157 

3  0 

5922.102643 

17.012997 

.462337 

.352620 

.654497 

E 

4  S 

230410.308065 

14.991034 

.461001 

.352992 

.646706 

S  0 

8533.044444 

16.766211 

.461002 

.352067 

.654226 

6  0 

5972.258448 

16.765768 

.464224 

.350912 

.658003 

7  0 

5880.106121 

17.001340 

.461848 

.352814 

.655003 

B  0 

5879.265444 

17.004706 

.461857 

.352801 

.654919 

9  0 

5879.264358 

17.004038 

.461852 

.352805 

.654923 

10  0 

5879.264165 

17.003685 

.461852 

.352810 

.654927 

• 

11  0 

5879.264128 

17.003843 

.461854 

.352812 

.654926 

12  0 

5879.264124 

17.003584 

.461856 

.352811 

.654928 

13  0 

5879.264124 

17.003645 

.461855 

.352811 

.65492B 

V 

IPAGE  4  BHDP  ESTIMATION  OF  CES  PRODUCTION  FUNCTION 

i 

THE  RESIDUAL  SUN  OF  SOUARES  ( 

*  5879.26 

)  NAS  SMALLEST  NITH  THE  FDLLDMING  PARAMETER  VALUES 

.* 

1  PI 

2  P2 

3  P3 

4  P4 

17.0036 

.461855 

.352811 

.654928 

a 

ESTIMATE  OF  ASYMPTOTIC  CORRELATION  NATRII 

PI  P2 

P3 

P4 

1 

2  3 

4 

PI 

1  1.0000 

» 

P2 

2  -.0021  1.0000 

P3 

3  .0973 

7663  1.0000 

P4 

4  -.9975 

0556  -.1251 

1.0000 

C 

THE  ESTIMATED  MEAN  SQUARE  ERROR  IS  534.5 

» 

►*- 

ESTIMATES  OF  ASYMPTOTIC  STANDARD  DEVIATIONS  OF  PARAMETER  ESTIMATES  NITH 

11  DEGREES  OF  FREEDOM  ARE 

\v 

1  PI 

2  P2 

3  P3 

4  P4 

3.27720 

2.604367E-02 

2.232392E-02 

3.2B3604E-02 

* 

IPAGE  S  BMDP  ESTIMATION  OF  CES  PRODUCTION  FUNCTION 

1 

D.6 

■; 

k  * 

>.\  V. 

*•».*  •/  \0  *„!  %" 

CASE 

OBSERVED 

PREDICTED 

COOK 

STD.  DEV. 

NO. 

RESIDUAL 

4  AF 

4  AF 

DISTANCE 

PREDICTED 

1  HI 

2  H2 

3  H3 

1 

-34.881964 

936.892380 

971.774344 

.355153 

12.727619 

326.0 

543.0 

753.0 

2 

29.828829 

883.645210 

853.816381 

.201149 

11.849887 

326.0 

543.0 

287.0 

3 

13.143130 

820.779310 

807.636180 

.040014 

11.934051 

326.0 

207.0 

753.0 

4 

-6. 158348 

669.880260 

676.038608 

.011607 

12.888390 

326.0 

207.0 

287.0 

5 

-6.039194 

837.946050 

843.985244 

.008278 

11.863760 

124.0 

543.0 

753.0 

6 

22.487536 

738.305280 

715.817744 

.148631 

12.751218 

124.0 

543.0 

2B7.0 

7 

.329927 

665.235910 

664.9059B3 

.000034 

12.919276 

124.0 

207.0 

753.0 

8 

-17.141728 

499.332690 

516.474418 

.062566 

11.642453 

124.0 

207.0 

287.0 

9 

17.808835 

893.147910 

875.339075 

.065858 

11.555498 

394.0 

375.0 

520.0 

10 

11.974304 

651.915660 

639.941356 

.041218 

12.675628 

56.0 

375.0 

520.0 

11 

-10.9082B1 

894.097370 

905.005651 

.025402 

11.651402 

225.0 

657.0 

520.0 

12 

.094793 

604.202090 

604.107297 

.000003 

12.640563 

225.0 

93.0 

520.0 

13 

-7.533425 

859.986050 

867.519475 

.011779 

11.553742 

225.0 

375.0 

911.0 

14 

-39.082710 

610.025140 

649.107850 

.441878 

12.697190 

225.0 

375.0 

129.0 

IS 

23.711366 

786.167230 

762.455864 

.020835 

6.060161 

225.0 

375.0 

520.0 

NUMBER  OF  INTE6ER  WORDS  OF  ST0RA6E  USED  IN  PRECEDIN6  PROBLEM  &71 
CPU  TINE  USED  .821  SECONDS 
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